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Abstract 

In these lectures I review observations of star-forming molecular clouds 
in our Galaxy and nearby galaxies to develop a physical intuition for 
understanding star formation in the local and high-redshift Universe. A 
lot of this material is drawn from early work in the field since much of 
the work was done two decades ago and this background is not generally 
available in the present literature. I also attempt to synthesise our well- 
developed understanding of star formation in low-redshift galaxies with 
constraints from theory and observations at high redshift to develop an 
intuitive model for the evolution of galaxy mass and luminosity functions in 
the early Universe. 

The overall goal of this contribution is to provide students with 
background helpful for analysis of far-infrared (FIR) observations from 
Herschel and millimetre/submillimetrc (mm/submm) imaging with ALMA 
(the Atacama Large Millimetre/submillimetre Array). These two 
instruments will revolutionise our understanding of the interstellar medium 
(ISM) and associated star formation and galaxy evolution, both locally and 
in the distant Universe. To facilitate interpreting the FIR spectra of Galactic 
star-forming regions and high-redshift sources, I develop a model for the dust 
heating and radiative transfer in order to elucidate the observed infrared (IR) 
emissions. I do this because I am not aware of a similar coherent discussion 
in the literature. 



8.1 Star- forming molecular clouds 

Here, I provide an overview of the observations and physics of star-forming 
molecular clouds in the Galaxy and nearby normal galaxies. The motivation 
is to develop the intuitive background for analysis of observations at high 
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redshift, hopefully avoiding naive and simplistic assumptions. The physics 

of the molecular ISM is not at all what one would have guessed before the 
advent of mm-line astronomy and the understanding developed over the last 
30 years is vital to galactic evolution studies. 

8.1.1 Background 

All known star formation in low-rcdshift galaxies occTirs in molecular gas 
clouds and inclusion of atomic Hi clouds in discussions of star formation 
is a 'red herring', except in as much as the atomic gas may feed the 
buildup of molecular clouds. In fact, there is no instance of star formation 
demonstrably observed to occur in atomic gas in nearby galaxies (as distinct 
from atomic gas dissociated by recently formed stars). Even in the extended 
ultraviolet (UV) disks of nearby galaxies, the young stars producing the 
UV could have been formed from molecular clouds simply too low in mass 
or surface density to be detected in present observations. The dominant 
process for formation of hydrogen molecules (II2) at low redshift is believed 
to be on the surface of dust grains where two H atoms remain in proximity 
long enough for a radiative decay to a bound state. The hydrogen molecule 
is dissociated from the ground state by a two-step process involving UV 
photon absorption in the Lyman and Werner bands longward of the Hi 
Lyman limit, followed by radiative decay to unbound vibrational states in 
the ground electronic state (HoUenbach et al. 1971). Thus, the destruction 
rate for H2 depends sensitively on the ambient UV radiation field and hence 
on the amount of dust extinction and H2 column density to the cloud surface 
- the latter since the foreground H2 may absorb all the incoming UV when 
the H2 Lyman and Werner bands become optically thick (i.e., self-shielding). 
Thus, whether a gas parcel is atomic or molecular will depend complexly on 
a number of factors: the volume density (since the formation rate depends on 
the densities of dust and atoms), the clumping of the ISM and the column 
density to the cloud surface and velocity dispersion which determine the 
dust continuum and H2 line opacities. For typical cloud densities of a few 
hundred per cm^^, the phase transition from Hi- to H2-dominated will occur 
at a column density A'^h+2H2 ~ 2 x 10^^ cm^^, corresponding to an extinction 
Ay = 1 mag for a standard dust to gas abundance ratio (Hollenbach et al. 
1971). 

The H2 molecule is homonuclear and therefore has no permanent dipole 

moment - its pure rotational transitions are weak quadrupole transitions 
with low emissivity. In addition, H2 has a small moment of inertia and its 
first rotational transition (J = 2-0) at 28.2 //m is therefore at relatively high 
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Fig. 8.1. The low-lying rotational levels of CO are shown on the left; the 
spontaneous decay rate for the J — 1-0 transition is 7 x 10~^ sec^^ and the critical 
density for pure coUisional excitation by H2 is ~3000cni~^. With line photon 
trapping and a typical optical depth of t ~10, the critical density for thermalisation 
is reduced to ~300cm~'^. Thus, the CO line brightness temperature provides a 
measure of the gas kinetic temperature. On the right, the CO and ^'^CO J = 1-0 
emission lines are shown in the direction of the Orion Kleinmann-Low (KL) nebula. 



energy {hv/k ~ 508 K), requiring a minimum gas temperature of Tk ~ 150 K 
for appreciable coUisional excitation. Thus trace molecules such as CO and 
HCN have become the standard surrogate probes of molecular clouds. CO 
and HCN have dipole moments of 0.1 and 2.3 Debye, respectively, and their 
larger moments of inertia (due to having heavier atoms than H) yield lowest 
rotational transitions at 115.2 and 88.6 GHz (2.6 and 3mm wavelength), 
corresponding to hu/k ^^5.5 and 4.2 K (see Fig. |8.1[ ). These molecules 
also have rare isotopic forms (e.g., ^^CO and C^^O) with slightly different 
moments of inertia, and hence shifted rotational lines which can be used to 
probe regions of high column density where the abundant isotope lines are 
very optically thick. 



8.1.2 Molecular excitation 

Excitation for the rotational levels of molecules like CO observed at mm and 
submm wavelengths is provided by collisions with H2. If these collisions are 
sufficiently frequent (compared to the spontaneous decay rate) the levels will 
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come into thermal equilibrium with the H2 and their excitation temperature 
will approach the gas kinetic temperature. For optically thin transitions, 
this occurs at the critical density nHj.crit = ^u-i/(<7fu-i) where ^u-i is 
the Einstein spontaneous decay rate (7 x lO^^sec"^ for CO, J = 1-0), and 
((Tt'u-i) is the collisional de-excitation rate averaged over the Maxwellian 
distribution of H2 thermal velocities. For CO (J = 1-0) this critical density 
is ~3000cm~^ and for HCN (and most other molecules) it is ~10^~^cm~^ 
due to their higher dipole moments and therefore higher A coefficients. 
Hence, the non-CO molecules are generally taken as probes of high-density 
molecular g compared with that traced by CO. 

Measurements of the rare ^^C isotopes of both CO and HCN indicate 
that the ^^CO emission is optically thick. The rare isotopes are typically 
seen at an intensity approximately 0.1-0.5 of the abundant species whereas 
the ISM abundance ratio is ^^C/C = 1/89-1/40, implying that the emission 
from the abundant species is optically thick and therefore saturated (see 



Fig. 8.1 right). For optically thick transitions, the upper-level population 
can be enhanced due to absorption of line photons, leading to excitation 
temperatures higher than those expected simply due to H2 collisions, since 
the line photons emitted upon spontaneous decay cannot easily escape the 
cloud. This so-called radiative trapping of the line photons builds up the 
radiation field at the frequency of the line, leading to enhanced excitation 
of the upper state via photon absorption. 

The escape probability formalism used to treat this optically thick 
situation was first applied to molecular clouds by Scoville &: Solomon 
(1974) and Goldreich & Kwan (1974) and is now routinely used to analyse 
interstellar molecule excitation (often called large velocity gradient [LVG], 
non-thermal equilibrium [non-LTE] analysis; see, e.g., van der Tak et al. 
(2007) and their pubhcly available RADEX code)[T} One of the biggest 
advantages of the LVG formalism is that it permits treatment of the coupled 
radiative transfer and molecular excitation as a local problem - a fact 
probably not fully appreciated by current routine users of these codes. This 
formalism is applicable to situations in which systematic velocity gradients 
are large compared to the small-scale thermal motions. The line photons 
from one region of the cloud are then incoherent with other regions due to 
the Doppler shift; they can then only interact with molecules in the local 
region near where they were emitted. 

In the photon trapping regime, the spontaneous decay rates {A) used 
in analysing the equilibrium molecular excitation are reduced by a factor /3 



t http : //www. strw.leidenuniv.nl/~moldata/radex.html 
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equal to the effective probability for escape of line photons from the emission 
region (Scoville & Solomon 1974; Goldreich & Kwan 1974). Thus, the 
critical density for thermalisation of the levels is reduced by the same factor 
p. For a spherically symmetric cloud with radial velocity field {ly (x r), the 
escape probability is given by 

/3 = -(l-e-0, (8.1) 
r 

where 

for T <^ 1, (3 = 1 and for r S> 1, /3 = l/r. Thus the critical density is reduced 
by a factor of r in the optically thick regime. For most giant molecular clouds 
(GMCs - see below), rco ^ 10 (based on the relative strength of ^^CO), so 
the critical density of CO (1-0) is reduced from 3000 to 300 cm~^. Thus, 
the CO excitation tem^perature will be approximately thermalised at the gas 
kinetic temperature (T\^) in most clouds. And since the CO line is optically 
thick, the brightness temperature at the emission line peak will indicate the 
H2 kinetic temperature. 

Similar considerations apply for the optically thick lines of molecules such 
as HCN and CS which have higher dipole moments - i.e., their critical 
densities are reduced by a factor l/r but since those critical densities are 
higher than that of CO, their excitation temperatures and line brightness 
temperatures will be subthermal (except in the densest cloud core regions). 
In fact, since the critical density for optically thick transitions is nH2,crit = 
i^u-i/ t) / {(^^u-i) J tlie radiative line strength (appearing in the 74-coefficient 
in the numerator and in r in the denominator) cancels out - i.e., the critical 
density then just depends inversely on the molecular abundance. 

In the regime of photon trapping, one should not assume that the emission 
of a rare isotope like ^^CO is optically thin simply because it is observed to 
have lower intensity than the abundant isotope emission. The rare isotope 
may simply have a lower excitation temperature since it has less photon 
trapping. This fact is under-appreciated by 'casual' observers. Typically, 
the ^"^CO emission is 1/5 to 1/10 of CO ~ does this mean the optical depth 
is always in the small range 1/5 to 1/10? - of course not! 

Molecular excitation can also be provided by absorption of continuum 
photons - either the IR continuum of dust or the cosmic microwave 
background radiation (CMB). In the escape probability formalism it is easy 
to include such continuum excitation. If (3 is the probability of a line photon 
escaping from the local region, then /3 is also the probability of an externally 
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produced background photon making it into the local volume. The low- lying 
molecular levels will maintain a base excitation temperature at the cosmic 
background level, 2.7 x (1+redshift) K (even in the absence of H2 collisions). 
However, since the line observations measure only the excess emission above 
the local continuum background, this CMB excitation does not result in 
detectable emission. The continuum from dust emission is not generally 
important for exciting low-lying molecular levels at mm wavelengths (since 
the continuum is optically thin and weak); but for higher-energy transitions 
at Arest < 500 fim the IR continuum should be included. 

8.1.2.1 Summary 
In summary, we can expect that: 

(a) The low CO levels will be thermalised and the brightness temperatures 
(as observed in spatially resolved clouds) will approximate the gas kinetic 
temperature, i.e., rB(CO) ~ Tk(H2). 

(b) Photon trapping in optically thick lines can greatly enhance the molecular 
excitation above the excitation due to collisions with H2 molecules. 

(c) The CO emission will trace the overall distribution of H2 for densities ~ 
10^""' cm^-^ ~ this turns out to be well tuned to characteristic densities of the 
majority of the molecular gas mass. 

(d) The emissions from higher dipole moment molecules (HCN, CS, etc.) will trace 
higher density regions with riHa ^ fO^~^ cm~^. 



8.1.3 Observed properties of molecular gas 

Extensive surveys of the CO and ^'^CO ( J=l-0 at A = 2.6 mm and J=2-l at 
A = 1.3 mm) emissions in the Galaxy have been used to determine the overall 
distribution of star-forming H2 gas and its properties (Scoville & Solomon 
1975; Sanders et al. 1985; Dame & Thaddeus 1985; Scoville et al. 1987; 
Clemens et al. 1988; Matsunaga et al. 2001; Dame et al. 2001; Jackson 
et al. 2006). CO emission is ubiquitous at all Galactic longitudes in the 
inner Galaxy (|/| < 90°) and within 1° latitude of the Galactic plane. The 
CO emission is particularly strong in the central 1° around the Galactic 
centre and within an annulus/ring from |/| = 20-50 degrees - the so-called 
molecular cloud ring at 3-7 kpc radius. This concentration of the H2 to 
the interior of the Galaxy is remarkably different from that of Hi; the 
21-cm Hi distribution is much smoother and the overall Hi distribution is 
fairly constant with Galactic radius out to i? ~ 12-15 kpc with mean column 
density perpendicular to the plane of (-/Vh) ^ 10^^ cm~^ (see Fig. 8.2). 
Within the Galaxy the overall mass contents of the molecular and 
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Fig. 8.2. The gas surface densities in the Milky Way disk are shown for H2 (Clemens 
et al. 1988), Hi (Burton & Gordon 1978) and giant Hii regions (>M42). Both the 
H2 and Hi values include a 1.36 correction factor for He. 

atomic gas are approximately equal (~ 2 x 10^ Mq in each) but most 
of the molecular gas is inside the Solar-circle whereas most of the 
Hi is outside 8.5 kpc. Since the total Galactic star formation rate is 
~3M0 per yr, the characteristic cycling time for H2 into new stars is 
~2 X lO^Mo/SMo/yr-lO^yr. Thus, the H2 clouds are forming stars at 
a rate much slower than their free-fall collapse time within the clouds 
((Gp)-V2^3-5Myr) but on a timescale much shorter than the Galactic 
lifetime. Star formation, is inefficient but a resupply of gas is required if star 
formation is not to die out in the next Gyr (which seems unlikely). 

The smaller-scale distribution of CO emission indicates that the H2 
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Fig. 8.3. The star-forming complex in Orion provides a close view of the complexity 
of both the ISM and the star formation activity. The left panel shows a very large- 
scale image of the FIR emission as imaged by the IR astronomical satellite (IRAS) 
(Courtesy NASA/JPL-Caltech). The upper-right panel shows the integrated CO 
(1-0) line emission from the two Orion GMCs (Ripple et al. 2012). A Hubble 
Heritage image of the visible Hll region M42 is shown at the lower right - the 
visible Hii region occupies a very small area of the Orion complex, as outlined in 
the IRAS image. 



gas resides in discrete clouds, rather than a diffuse, continuous medium. 
Along a typical line of sight in the Galactic plane there will be three 
to six kinematically discrete CO emission features, occupying a small 
fraction of the overall permitted velocity range. Spatial mapping of the 
individual emission features indicates extends perpendicular to the line 
of sight typically 10 to 80 pc - hence the name. Giant Molecular Clouds 
(GMCs). 

The Orion CMC (actually two clouds) is the nearest CMC with an overall 
extent of 40 pc (see Fig. 8.3) - in striking contrast to the optically visible 
Orion Nebula (M42) for which the extent is ^Ipc! In the vicinity of M42, 



the CO brightness temperature increases to a peak of ~60K (see Fig. 8.1 



right), compared to 5-20 K in most of the rest of the cloud. These elevated 
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Fig. 8.4. The internal velocity dispersions of GMCs arc shown as a function 
of diameter for clouds with and without giant Hii regions, i.e., Hll regions 
more luminous than M42 (Scoville et al. 1987). This illustrates the empirical 
correlation found between cloud size and lincwidth (the so-called sizc-lincwidth 
correlation). The Hii region clouds depart from this size-linewidth at low 
masses/sizes presumably due to feedback effects from massive star formation. 

temperatures are largely due to heating by the dust-embedded cluster of 
young stars in the Kleinmann-Low (KL) IR source, located behind the visible 
Trapezium/M 42 OB star cluster. 

The gas heating associated with the KL nebula (and other IR sources 
associated with active star formation) occurs in two steps: heating of the 
embedding dust by absorption of photons emitted by the young stars, 
followed by H2 heating through collisions with the heated dust grains. 
Thermal equilibrium between the dust and the H2 is theoretically expected 
to occur at gas densities riHa > lO^cm"^ (Goldreich Sz. Kwan 1974). Indeed, 
the dust colour temperature in KL is ~60K, i.e., similar to the gas kinetic 
temperature determined from the CO line brightness temperature. 



One of the most perplexing and still not understood features of the 
GMCs is their large velocity dispersions, as determined from the 
width of the molecular emission lines. At gas temperatures of 10- 
80 K, the thermal velocity dispersion (sound speed) should be ~0.1- 
0.3km/s, yet the observed velocity dispersions are typically l-4km/s. 
The linewidths are correlated with cloud diameter, D. The observed 
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relation is ay = 0.31Dp^^^ km/s (see Fig. 8.4). Thus the kinetic 
energy in large-scale supersonic motions is ~100 times the expected 
thermal energy. The containment of these supersonic turbulent motions 
is provided by the self-gravity of the molecular gas within each cloud. 
Independent estimates of the cloud masses using both dust extinction 
measurements and non-LTE analysis of the molecular lines, indicate 
approximately enough gas mass for the GMCs to be self-gravitating and 
in virial equilibrium with the observed internal motions. Nevertheless, it 
should be pointed out that although the kinetic energy and gravitational 
potential energy approximately balance, the GMCs do not have the expected 
spherical shapes; they are often elongated and have internal filamentary or 
sheet-like structures. 

In any case, it is clear that the GMCs are not in pressure equilibrium 
with the external diffuse Hi and Hii phases of the ISM. Their pressures are 
generally taken to be riT = 3000 cm~^ K, i.e., a factor of 100 lower than the 
turbulence pressures within the GMCs as reflected in their internal velocity 
dispersions. 

Given that the GMCs are self-gravitating, one can obtain their mean 
densities as a function of cloud diameter using the above empirical relation 
between the velocity dispersion and cloud diameters (i.e., the observed size- 
linewidth relation): 

(nn,) = 180 i^—j cm-3, (8.3) 

where the scaling is to 40 pc - the size at which half of the overall Galactic 
H2 mass is in clouds larger and half is in clouds smaller. The cloud mass 
distribution function is N(M) (x M^^'^ . For a cloud of diameter 40 pc, 
the above relation implies a mass of 4 x 10^ Mq (including the 36% mass 
contribution of He) . 

8.1.3.2 Molecular masses from r ^ 1 CO emission 

How is it that the emission in optically thick CO lines can be used to estimate 
the total mass of H2 in galaxies? Observers of extragalactic CO commonly 
make use of a so-called X-factor or aco to translate measurements of CO 
emission into estimates of molecular gas masses - let us see how this is 
physically justifiable. 

For a resolved cloud the integrated line brightness temperature, Leo = 
/ Tedf;, can be integrated over the projected area of the cloud to yield a 
'luminosity' Leo = d'^ J Icod^ where d is the distance to the cloud. Thus, 

Leo = Tb{C0)AV7tR^ K km/s pc^, (8.4) 




Fig. 8.5. The CO luminosities of GMCs are shown as a function of their virial 
masses for clouds with and without giant Hii regions (Scoville et al. 1987). 



where R is the radius of an assumed spherical cloud. For clouds in virial 
equilibrium, Ay = (GM / RY/"^ and therefore, 

Leo = (37rG/4p)V2TB(CO)MGMC, (8.5) 



where p is the mass density. Equation |8.5| predicts a linear scaling between 
cloud mass and the CO luminosity provided the clouds have approximately 



similar mean density and temperature. Equation 8.5 indicates that 
the constant of proportionality between Lqq and M will vary as T j ^fp. 
Physically, what is going on here? Although the CO line is optically thick 
(and often thermalised) , the increased mass of larger clouds gets reflected in 
an increased surface area emitting CO photons and an increased linewidth 
over which they are emitted. 



Figure 8.5 illustrates the extremely good empirical correlation between 
measured CO luminosities and their virial masses (determined from resolved 
measurements of the size and linewidth). In fact, if the temperatures are 
determined by radiative heating of dust within the cloud (plus a low-level 
background heating by cosmic rays heating the gas to ~ 10-15 K), the overall 
range of mean temperatures is expected to be rather small. The dust 
radiatively cools as tx and it is hard to get the dust much above 

40 K except in very localised regions of active star formation. Typical FIR 
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colour temperatures of GMCs are in the range 20-40 K and only a few ultra- 
luminous IR galaxy (ULIRG) nuclei get up to 60 K. Regions with elevated 
gas temperatures and hence high rB(CO) are also likely to be denser (since 
their high temperatures probably reflect more concentrated star formation 
due to higher gas density) . The effects of increased temperature and density 



will therefore partially compensate each other (see Equation 8.5). 



For extragalactic observations where individual GMCs are not resolved, 



Equation 8.5 can still be used since the total CO luminosity is just the sum 
of that of the individual clouds, provided the GMCs do not overlap both on 
the sky and in velocity (unlikely given their small-volume filling factor in 
normal galaxies). In the extragalactic observations, the cloud diameters are 
not usually measurable but the unresolved apparent brightness temperatures 
vary as distance d"^ so that the distance to the galaxy may replace R in 



Equation |8.4[ For the external galaxies, Equation |8.5| is translated to 

Mu,{Mq) = acoLco K km/s pc^, (8.6) 



where aco — 4.9 (Solomon & Barrett 1991). 

We should also note that for modest changes in metallicity (affecting the 
CO abundance relative to H), qco will change slowly (until CO becomes 
optically thin in large areas) - this is dramatically demonstrated by the 
fact that ^^CO luminosities of clouds are typically only a factor of 4-10 
times lower than those of CO despite the much lower abundance ratio of 
^^C/C (1/40— 7>l/90). For the optically thick regime where photon trapping 
is important for the CO excitation, one can show analytically that the aco 
should scale as CO abundance or metallicity Z"'''^ (see Scoville & Solomon 
1974). 

One instance where the above analysis must be modified is in a galactic 
nucleus if the gas is smoothly distributed in the central galactic potential. 
Here, the CO linewidth will be determined not by the self-gravity of indi- 
vidual clouds but instead by that of the stars plus the gas, i.e., the linewidth 
will be larger than that associated with just the gas mass. There will then 
be more CO photons emitted per unit gas mass. For parameters associated 
with the most extreme cases in ULIRGs, the CO conversion factor may be 
reduced by factors of 2-5 (see Downes & Solomon 1998 and Bryant & Scov- 
ille 1999). In the photon trapping regime, the CO excitation temperature 
(and hence the brightness temperature) of the line will vary as the molecular 
abundance to the ~0.4 power (Scoville &: Solomon 1974). Thus, one expects 
that the CO line emissivity to mass conversion factor will vary with metal- 
licity as Z~^''^ at high redshift where the metallicities are generally lower. 
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8.1.3.3 Lifetimes of GMCs 

Although much current theoretical analysis of star formation adopts the 
point of view that GMCs are relatively short-lived (10-30 Myr), I would like 
to provide several arguments why this can't be the case. These arguments 
were well known in the early days of molecular line astronomy but seem to 
have been overlooked in current discussions. 

The first argument is simply based on the required conservation or con- 
tinuity of ISM mass between the various ISM phases (Hi, Hii and H2). If 
the GMCs, or more specifically the H2 molecules, are short-lived then there 
must be cycling of gas between these phases. The mass flux from one phase 
to another and back must then be in equality (see Scoville & Hersh 1979; 
Koda et al. 2009). This mass flux equality is given simply by 

A^H2^-HI+HII (= MYi^/Tn^) = -^HI+HII-^H2 (= ^HI+HIi/thI+HIi) , (8.7) 

where we assume relatively little cycling to young stars (since that timescale 
was found to be ~1 Gyr above). In the interior of our Galaxy and most 
nearby spirals, the H2 dominates in mass compared to Hi and Hil by a 
substantial factor (e.g., more than a factor of five in the Galactic molecular 
cloud ring) . Thus if cycling between the phases is to be plausible at all (with- 
out disobeying mass conservation!) at these interior radii, the timescale for 
a molecule to remain a molecule must be correspondingly longer than the 
timescale for a particle in the diffuse phases (Hi or Hii) to remain atomic or 
ionic. The latter might be expected to be given by the dynamical time at 
density a few H cm~^ (~10^ yr) or by the time needed to pass between the 
compressive spiral arms (also ~2 x 10^ yr). Inverting the mass continuity 



equation (Equation 8.7), one finds: 



= Jf^' X THI+HII ~ 5 X (2 X 10^) ~ 10^ yr. (8.8) 
jWhi+hi 

Thus, for cycling between the phases to work, the characteristic lifetime of 
a H2 molecule should be at least several times 10^ yr. And if there is no 
such cycling, then the timescale is even longer. Although this calculation 
was done for the interior of the Galaxy where H2 dominates Hi, the estimate 
can be expected to hold elsewhere since it is in the interior of the Galaxy 
where there is the most disruptive feedback from star formation and other 
processes. In addition, the internal structure of the GMCs in the outer 
Galaxy is probably similar to that in the inner Galaxy (i.e., they are equally 
difficult to disrupt). 

One must note that this is the timescale for a molecule to remain a 
molecule and not necessarily the timescale for a given CMC to retain its 
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Fig. 8.6. The scaleheights of the GMC distribution perpendicular to the Galactic 
disk are shown as a function of cloud virial mass (Scoville et al. 1987). The nearly 
square-root dependence on cloud mass suggests an equipartition of the cloud kinetic 
energy distribution, requiring that the GMCs last longer than several 10* yr. 



physical structure. Clouds may and probably do grow within the spiral arms 
and fragment into smaller clouds outside the arms - but the molecules cannot 
cycle all the way back to the diffuse Hi or Hii phases on a short timescale. 

Another line of reasoning suggesting long GMC lifetimes is the fact that, in 
general, feedback processes from star formation within the GMCs are ineffi- 
cient at disrupting the massive clouds (Scoville & Hersh 1979; Scoville 2003; 
Murray et al. 2010). Most of the stellar feedback will be radiated and dissi- 
pated in shock fronts or leak out the cloud boundaries (as is surely the case 
in M42). The total momentum needed to unbind the clouds is enormous, 
typically ~ 5-10 x lO^MQkm/s. Lastly, I point out that the scaleheights 
(perpendicular to the Galactic disk) of the GMCs as a function of their mass 
suggest that they have achieved approximate equipartition. In Fig. |8.6[ the 
2;-scaleheights of clouds indicate that oc M^^''^, i.e., approximately the 
1/2 power expected for equipartition. Equipartition requires that the clouds 
survive at least several cloud-cloud scattering times. The latter varies with 
the local Galactic space density of clouds but is almost always > 10^ yr. 
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8.1.3.4 GMC supersonic internal motions 

I mentioned earlier that the GMCs have highly supersonic internal motions 
as determined from the molecular line Doppler widths and that to date there 
is no satisfactory explanation. The long cloud lifetimes highlight this prob- 
lem. The maps of molecular emission within individual clouds show that 
the molecular emissions and their large velocity spread are fairly smooth 
across the projected area of the clouds, implying that the area covering fac- 
tor is > 1. If this is the case then the supersonic turbulent gas parcels will 
necessarily collide and dissipate their energy with the time it takes to cross 
the cloud diameter. This timescale is 40pc/3km/s ^■lO'^yr. Thus the tur- 
bulent energy required to support the GMCs against gravitational collapse 
needs to be replenished within a similar timescale. (Reduction of the dissi- 
pation by postulating that magnetic fields are driving these motions requires 
very strong fields and does not change this estimate significantly since the 
collisional dissipation can still occur along the field direction.) 

Although the feedback energy released from embedded young stars is of 
the right order of magnitude to replenish the turbulent energy, most of 
the energy released from the protostcllar outflows is radiated away in high- 
velocity shock fronts and is also deposited on small length scales rather than 
scales comparable to the cloud size needed for the large scale support. Tur- 
bulent cascades generally transport energy from large scale to small scale 
(except in 2D turbulence, see Robertson &; Goldreich 2012) so it is difficult 
to see how the stellar feedback can maintain the turbulent energy. 

A possible source of the turbulent energy might be large-scale 'corruga- 
tions' in the galactic plane mass distribution. It should be pointed out 
that despite the apparent virial equilibrium of the GMCs (in the sense that 
their gravitational potential energy is approximately twice their internal ki- 
netic energy), the morphologies of the clouds do not appear virialised and 
spherical - this provides a strong argument that they are constantly being 
distorted by the external force gradients. Those associated with normal den- 
sity wave spiral structure pass by only every ~ 2 x 10^ yr and are therefore 
not sufficient. Jog & Ostriker (1988) propose that GMC-GMC scattering is 
viable to maintain the internal motions but such scattering occurs also on a 
typical timescale of only every ~ 10^ yr, which seems too long compared to 
the dissipation timescale. 

D'Onghia et al. (2011) have recently suggested that many observed spiral 
structures may be stochastic in nature, where mass seeds corresponding to 
GMC masses (very likely the GMCs themselves) will induce a mass overden- 
sity in the stellar disk via the swing-amplification process. Is it viable that 



16 



Nick Z. Scoville 



these associated stellar disk mass enhancements actually feed back gravita- 
tional perturbation energy into the internal supersonic motions of the GMC 
gas? 

8.1.3.5 Summary 

In summary, the observations of Galactic GMCs and local galaxies indicate 
that : 

(a) The CO emission arises from discrete, self-gravitating clouds (GMCs). 

(b) These GMCs have high internal supersonic motions such that their overall ki- 
netic energy content is dominated by turbulence with effective pressure 100 
times the gas thermal pressure and 100 times the pressure of the external 
diffuse ISM. 

(c) The GMCs have a mass spectrum N{M) oc M~^-^ with the midpoint in the 
mass contributions at cloud diameter ~ 40 pc and mass 4 x 10^ Mq. For this 
size cloud, the mean density is riHs ~ 180 cm~^ - larger clouds have lower 
density and smaller clouds higher mean density. 

(d) For such self-gravitating clouds with optically thick CO Uncs, the emitted CO 
luminosity is approximately proportional to the cloud mass with a constant of 
proportionality scaling as T/(p)^/^ (i.e., aco varying as p^^"^ /T) and having 
much lower than linear dependence on metallicity. 

(e) The hfetimes of the H2 and probably the GMCs are apparently larger than 
10® yr and possibly 10^ yr based on simple mass-conservation arguments in the 
ISM between the diffuse and dense phases. This is consistent with what one 
might expect based on the great inertia of the GMCs (i.e., resistance to disrup- 
tion/dissociation) and the fact that their extraordinarily high effective internal 
pressures (compared to the external diffuse phases) makes it very difficult for 
pressme disturbances in the external medium to significantly influence their 
internal structure. 

(f) The maintenance of the supersonic motions within GMCs, providing support 
against collapse, remains unsolved since this energy must be replenished on the 
timescale of a few 10^ yr, yet the clouds last > 10® yr. Internal feedback from 
star formation deposits energy on length scales which are too small. The non- 
spherical shapes of the GMCs suggest large-scale external force gradients may 
be responsible, such as those associated with GMC-GMC scattering, clouds 
motions perpendicular to the galactic disk potential or the recently proposed 
stochastic spiral structure. 



8.2 Star formation 

Here I describe the various probes of star formation in galaxies and their 
correlation with the molecular gas contents. I devote considerable effort to 
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developing an interpretive framework for the FIR emission from optically 
thick dust clouds. 



8.2.1 Probes of star formation 

There are a number of observational tracers which have been developed to 
measure the star formation rates (SFRs) in galaxies (see Calzetti, this vol- 
ume). The Hi recombination lines in the visible and near-IR (e.g., Ha and 
Pa) have the fluxes proportional to the Hii region emission measures, hence 
the OB star formation rate over the last 10'' yr. The restframe (far-)UV 
continuum ([F]UV), at A < 2000 A, arising from hot, early-type stars has 
been used to infer the SFRs for large samples of galaxies. Both the emission 
lines and the UV continuum can be severely attenuated by dust extinction 
in star-forming regions. Even for the galaxies with detected UV continuum, 
the extinction corrections are often factors of 5-10! For the dust-obscured 
star formation, the FIR luminosity (Arest = 8 — 1000 fim) provides a much 
more reliable measure of the SFR. These FIR SFRs can now be obtained for 
large samples of galaxies using observations from the Spitzer and Herschel 
space telescopes, albeit with relatively low angular resolution and sensitivity 
to SFR (compared to the UV). A summary of all these techniques, including 
relevant SFR equations, appears in the recent paper by Murphy et al. (2011) 
so I will not detail all of them here - instead I will focus on developing a 
physical understanding of the IR emission. 



8.2.2 Infrared emission 

The FIR emission from both star-forming and active galactic nucleus (AGN) 
sources arises from dust surrounding these sources which has been radia- 
tively heated by absorption of the outflowing photons. Here I develop the 
logical steps for interpreting the IR emission since I have not seen this done 
systematically elsewhere. 

The dust temperatures are determined by radiative equilibrium at dis- 
tance R from a central source of luminosity L, with 

47r4(e^)ard^ = ^al{K^)L/ {A-kR^), (8.9) 

where Od and are the dust grain radius and temperature, and {e^) and 



{Ky) are the dust emission and absorption efficiencies (shown in Fig. 8.7 left 
panel), weighted, respectively, by the Planck spectrum at the local dust tem- 
perature and that of the luminosity source heating the dust. If the emission 
efficiency varies as 1/A, i.e., ej^ oc T^, then 




Fig. 8.7. Left panel: the dust absorption opacity as a function of wavelength (Isella 
et al. 2010) for standard interstellar grain composition (12% silicates, 27% organics, 
and 61% water ice; Pollack et al. 1994), and size distribution n{a) cx a""^-^ 
from 0.01 to 1 /im radius. The absorption coefficient is normalised relative to 
that at visual wavelength (6060 A, ny — 2.3 x lO^gr"^) - for a standard gas- 
to-dust abundance ('^lOO in mass) with A'^h-i-2H2 = 2 x 10^^ cm~^ per mag of visual 
extinction {Ay). An important feature of the extinction curve is the relatively 
flat broadband absorption coefficient over the range A = 2-50 /im. In the right 
panel, I show the Planck integrated absorption/emission coefficient as a function 
of blackbody temperature. Due to the 'flat' absorption coefficient in the near- 
and mid-IR (NIR/MIR), the Planck-integrated absorption and emission efficiency 
is quite independent of temperature from 1000 K down to 100 K. The dust opacity 
curve was calculated by Andrea Isella and is available from him or myself. 



LV5 

Td oc ^ (8.10) 



(Goldreich & Kwan 1974). More generally, 



I (£„(Td)) 



where {Ky{T)) and {e,y{T)) are shown in the right panel of Fig. 8.7, Due 
to the flatness of the broadband absorption coefficient from 2 to 50 /xm, the 
grain absorption and emission efficiencies integrated over black-body spectra 



are decreasing only modestly from Tbb =1000-100 K (Fig. 8.7, right). 

For an optically thin dust distribution surrounding a luminous source, the 
dust heating is mainly due to the central short-wavelength source, but for 
an optically thick dust envelope, the interior dust does not see the central 
source and is instead heated by secondary radiation. This secondary radia- 
tion, having longer wavelength than the central stellar or AGN source, is less 
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efficiently absorbed and the dust in the opticahy thick case wih therefore 
be colder (than it would be if exposed to the shorter-wavelength primary 
photons of the central source). In very optically-thick cases, the {ki,{Tl)) in 
Equation 8.8 can be evaluated approximately with Tl{R) ~ T^{R), appro- 
priately weighted over nearby radii. 

Figure |8.8| shows the computed dust temperatures for the optically thin 



and optically thick cases, evaluated from Equation 8.11 In the optically- 
thick regime, the fall-off in dust temperature is oc r~^/^ - an extremely 
simple form due to the fact that the photons heating the grains and those 
emitted by the grains have similar wavelength distributions. One will there- 
fore have (k) ~ (e) in the very optically-thick dust clouds. 

The radial scale in Fig. 8.8 is for a central source luminosity of 10^^ L©, 
appropriate to ULIRGs and submm galaxies (SMGs), but the radial dis- 
tances can be scaled as L^/^ for other luminosities (see Equation 8.11). 
Thus, these temperature profiles can be equally well applied to a dust cloud 
surrounding a luminous protostellar cluster of luminosity 10^-10® Lq. The 
temperature profiles in Fig. |8.8| start at ~1000K which is a little below the 
dust sublimation temperatures (< 1500 K). Inside this radius, the dust will 
not survive. For a less luminous source, this inner radius will scale inwards; 
but the physical scales will all change by the same L^/^ and the modelling 
remains homologous. 



8.2.3 Dust optical depth: r<l or t > I? 

It is quite trivial to observationally distinguish the optically thick and thin 
sources since the latter will have a power-law flux distribution on the short- 
wavelength side of the peak whereas the former will appear more exponential 
(see Scoville & Kwan 1976). (If the dust distribution is clumpy, photons can 
emerge from the inner regions with hot dust - thus, in optically-thick dust 
clouds, a non-exponential spectral energy distribution [SED] at short wave- 



lengths might also occur; see Fig. 8.17 ) 



Virtually all IR sources associated with active star-forming regions (Galac- 
tic GMCs and starburst nuclei) are optically thick into the MIR (based on 
their sharp short wavelength fall-off). Since the dust opacity is fairly flat 
across the wavelength range (2-50 /^m), one then expects that the source 
will become optically thin only at A > 50 /im, as long as it is optically thick 
at 3-10 /im. Thus, it is most appropriate to employ the optically thick dust 



temperature distribution shown as a solid line in Fig. 8.8 and fit numerically 
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Fig. 8.8. The temperature of dust heated by a central luminosity source is shown 
as a function of radius for optically thin dust (with all the heating provided by a 
~ 10* K blackbody) and for very optically thick dust (where after the innermost 
radius, the heating is due to dust at the same radius and temperature). This 
resembles closely the very optically-thick case since the temperature gradients are 
generally quite shallow. It is important to appreciate that the radial physical 
scalelengths will simply stretch homologously with changing source luminosity. 
Thus the same curves can be used for sources with very much higher or lower 
luminosity. For optically-thick FIR dust emission the dust temperature varies 
as Td oc r^^/^. In real emission sources, the dust is optically thin right at the 
inner boundary, optically thick at intermediate radii and then thin at the outer 
radii so the temperature profiles must be pieced together at the appropriate radii. 
(These equilibrium dust t emp eratures were calculated using the grain absorption 
coefficients shown in Fig. |8.7| and the inner boundary is taken to be where the 
grain temperature is ~ 1000 K, i.e., somewhat below the expected dust sublimation 
temperatures.) 
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by 

Td = 630^^ -j:^ K. (8.12) 



One can invert Equation 8.12 to find the characteristic size of the emitting 
region: 

i?thick = 100 (63/rd)' (L/IQI^Lq)'/' PC (8.13) 

if the total IR luminosity and dust temperature are known. The latter might 
be derived by fitting the MIR SED, or more crudely, from the wavelength of 
the peak. If the dust is opaque, then the standard Planck expression yields 
Td = 100 X (51 ^m/Apeak) K but if the dust is optically thin, then the peak 
wavelength is reduced by a factor ~ (3/(3 + a)) where a is the power-law 
index for the opacity, oc i^", near the IR peak. 

In practice, the opacity can never be much greater than unity at the peak. 
This is due to the simple fact that the luminosity cannot escape from the 
inner regions which have high opacity to the cloud surface. And once the 
outward luminosity flux has shifted to wavelengths where the opacity be- 
comes less than unity, it escapes. This can be seen analytically using the fact 
that the emergent emission for a given grain is proportional to Bu€ue^~'^''\ 
At A > 70 ^m, Ey ccu"""^-^ . If r^/ = (z^/z^o)^'^ (i-e., unity at v = vq), then the 
emission of the grain will peak at z^m — (3 + 1.6)z^o/(l-6 + huo/hTf^), where 
T = {i^m/i'o)^'^ ■ Thus, for example, with Td = 100 K, the blackbody peak at 
51 //m is shifted to 78 and 395 ^m for Aq = 102 and 787 ^m, respectively. At 
the peaks the optical depths are Tpeak = 1-5 and 3.0 in the two cases. This 
illustration was for an isothermal dust distribution with r simplistically rep- 
resenting the foreground optical depth. This also provides a strong note of 
caution - if one is deriving the dust temperature from the wavelength of 
peak emission, the peak will usually be shifted to a longer wavelength where 
the dust starts to become transparent. 



8.2.4 Dust temperature of the emergent luminosity 

So how can the characteristic dust temperatures be derived? Probably the 
best approach, if the dust is believed to be optically thick on the short wave- 
length side of the peak (based on an 'exponential' rise), is to simply assume 
that r ~ 1 at the peak and estimate the temperature of the dust emitting 
the bulk of the emergent luminosity from Td > 100 x (80/im/Apeak) K- To 
improve on this ad hoc correction requires an assumption of the dust density 
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Fig. 8.9. The conversion of observed IR luminosities into estimates of the SFR 
depends on the duration time of the dust absorption. Here I have used Starburst99 
models with a Kroupa IMF and assumed that all stellar and nebular photons 
longward of the Lyman limit are absorbed for the dust duration time (in Myr) and 
then none are absorbed after that. For normal star-forming regions this timescale 
is --10 Myr and for ULIRGs 50-100 Myr. If the dust envelope lasts less than 10 Myr 
and or is only partially covering, then the conversion factor is significantly higher 
- both are certainly true for an exposed Hii region like the Orion nebula. 



distribution with radius and modelling the emergent SED (as is done below 
in Section 8.2.9). 

For an optically-thick source, the emergent radiation at each wavelength 
will be from a depth where tx — 1- Since the dust opacity falls off at longer 
wavelengths and the temperature is falling at larger radii, this implies, 
somewhat counter-intuitively, that at wavelengths short of the peak wave- 
length, one will sample dust at increasingly larger radii, i.e., lower and lower 
temperatures! 



8.2.5 Star formation rate from Ljr 

Derivation of SFRs from the IR luminosities is fairly straightforward and ro- 
bust. This was first done by Scoville & Young (1983) using the fact that the 
bulk of the luminosity from a stellar population at early times is generated 
largely by the OBA stars. For those stars an approximately fixed percentage 
(13%) of their initial mass gets processed through the CNO cycle and one 
knows the energy produced per CNO process. This derivation is analogous 
to the 'fuel consumption theorem' of Renzini & Buzzoni (1986). 

A more precise, modern approach is to run starburst models (e.g., Star- 
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burst99; Leitherer et al. 1999) and integrate up the luminosity as a function 
of time. One must assume a stellar initial mass function (IMF), and its 
mass range, and decide which photons will be absorbed by dust and for 



how long the dust absorption persists. In Fig. 8^ I show results obtained 
using Starburst99 with a Kroupa IMF (0.1 to 100 M©, Kroupa 2001) as a 
function of the duration of the dust absorption. The latter quantity is prob- 
ably ~ 10'' yr for Galactic star-forming regions but could be closer to 10^ yr 
for merging starburst galaxies where the dust is more widely distributed. For 
this calculation, I assume all stellar and nebular photons longward of the 
Lyman limit (912 A) are absorbed by the surrounding dust. From Fig. 8.9 
one could reasonably compute a SFR given by 

SFRiR = 2 - 1.2 X 10-1° {Lik/Lq) Mq/jv, (8.14) 

with the lower value being appropriate to the ULIRGs which have longer 
duration for the dust shrouding. The standard relation given by Murphy 
et al (2011) corresponds to 1.5 x IO^^^Mq yr^^ (Ljr/Lq). The simpler 
derivation outlined above (based on the CNO cycle energy production) is 
quite similar to Equation |8.14| after one corrects for the mass going into 
non-OBA stars for a Kroupa IMF. 



8.2.6 Dust and ISM mass estimates 

On the long-wavelength Rayleigh-Jeans (R-J) tail of the FIR emission, the 
dust will be optically thin and the observed continuum fluxes provide an 
excellent means of determining the overall mass of dust. If the dust-to-gas 
abundance is normal, this dust mass can then be scaled to estimate the 
overall mass of ISM within a star-forming region or a distant galaxy. 

On the optically-thin R-J tail of the IR emission, the observed flux density 
is given by 

= K,rdusti^'Mdust/(47rdf), (8.15) 
or in terms of the dust opacity per unit gas mass, Ku{ISM) = Ki, x 

Mdust/MisM, 

= K,(ISM) rd,sti^'MisM/(47r(if), (8.16) 

where di is the source luminosity distance. In normal star-forming galax- 
ies, the majority of the dust is at ~20-25 K, and even in the most vigorous 
starbursts like Arp220 the FIR/submm emission is dominated by dust at 
temperatures < 45 K. Thus the expected variations in Tdust have less than a 
factor two effect on the observed flux. 
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Fig. 8.10. Left panel shows the dust-to-gas mass ratios derived by Draine et al. 
(2007) for galaxies from the SINGS nearby galaxy survey, selecting only those 
galaxies with both SCUBA 850 /im fluxes and complete maps of the H2 and 
Hi gas; right panel shows the mass-metallicity relation for low-z galaxies (grey 
points) and binned values for z~2 galaxies (red dots) from Erb et al. (2006). 
Over a range of 0.5 dex in metallicity below that of the Milky Way there is little 
variation in the dust-to-gas ratios. Since the galaxies selected here are massive, their 
metallicities even at z = 2 are expected to be within this range based on emission 
line ratios. Lower-metallicity irregular galaxies probably do show a decrease in the 
dust abundance (see Draine et al. 2007). 



The dust opacity per unit mass of total ISM gas, Ki,(ISM) in Equa- 
tion 8.16 , can be calibrated from the extensive submm observations of nearby 
galaxies. Seventeen of the nearby SINGS survey {Spitzer Infrared Nearby 
Galaxies Survey; Kennicutt et al. 2003) galaxies have good total submm flux 
measurements obtained with the SCUBA instrument, mounted at the James 
Clerk Maxwell Telescope, at 850 /xm (see Draine et al. 2007), as well as good 
measurements of the total molecular (H2) and atomic (Hi) gas masses. 

In Fig. |8.10 the derived dust-to-gas (H2 + Hi) mass ratios from Draine et 
al. (2007) are shown for the galaxies having SCUBA 850 fim measurements, 
for a range of spiral type (Sa to Sd) and as a function of mean metallicity. 
(Equivalent data for low-redshift elliptical galaxies are not available from 
Draine et al. 2007 due to lack of gas mass measurements in the early type 
galaxies.) 

Figure 8.10 shows that over a range of ~0.5dex in metallicity, there is lit- 
tle empirical evidence of variation in the dust-to-gas mass ratios and hence 
the submm flux. If one includes even lower-metallicity galaxies that do not 
have SCUBA 850 fim fluxes (hence the submm fluxes must be extrapolated 
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Table 


3.1. Low-z 


galaxies 


with submm 


& ISM data 


Galaxy 


Distance 


S,(450Ai) 






Ids' Mh 




(Mpc) 


(Jy) 


(Jy) 


{Me) 


(Mo) 


NGC4631 


9.0 


30.7 


5.73 


9.2 


9.5 


NGC 7331 


15.7 


18.5 


2.98 


9.4 


9.7 


NGC 7552 


22.3 


20.6 


2.11 


9.7 


10.0 


NGC 598 


76.0 


2.3 


0.26 


9.8 


10.1 


NGC1614 


62.0 


1.0 


0.14 


9.7 


10.0 


NGC 1667 


59.0 


1.2 


0.16 


9.3 


9.6 


Arp 148 


143.0 


0.6 


0.09 


9.9 


10.2 


1ZW107 


170.0 


0.4 


0.06 


10.0 


10.3 


Arp 220 


79.0 


6.3 


0.83 


10.0 


10.3 


12112+0305 


293.0 


0.5 


0.05 


10.3 


10.6 


Mrk231 


174.0 


0.5 


0.10 


9.8 


10.1 


Mrk273 


153.0 


0.7 


0.08 


9.9 


10.2 



from an overall SED fit using shorter- wavelength observations) , there is evi- 
dence for a metallicity dependence in the dust-to-gas ratio. Lastly, it should 
be emphasised that although some variations in the ratio of submm flux 
to ISM mass may be expected, ISM mass estimates at ~30% accuracy (see 
Fig. 8.10) are still likely to be at least as accurate as those from CO line 
measures and much quicker, enabling large samples to be analysed. 

In view of the large uncertainty in the submm absorption coefficient k and 
its scaling with frequency, we adopt an empirical approach based on submm 
observations of local galaxies where Hi and H2 masses have been estimated. 
For the local galaxies on which to base this empirical approach, it is vital 
that both the submm fluxes and ISM masses are global values. In addition 
as a check on the reliability of the submm measurements we require two 
long wavelength flux measurements so one can check if there is reasonable 
consistency with expected values of the spectral index /3. 

For the spectral index /3 of the R-J tail {S^ varying as u^), the observed 
flux ratios of submm galaxies can vary between 3 and 4. For most dust mod- 
els the spectral index of the opacity is typically 1.5 to 2, implying /3 = 3.5 to 4. 
Empirical fits to the observed long wavelength SEDs give /3 = 3.5 to 4 (Dunne 
& Eales 2001; Clements et al. 2010) for local galaxies. For high- 2; submm 
galaxies, the spectral index can be between 3.2 and 3.8, but in most cases 
the shorter-wavelength point is getting close the IR peak in the restframe 
and therefore not strictly on the R-J tail. In the following, we adopt /3 =3.8. 



Table 8.1 lists local spiral or star-forming galaxies for which both 450 and 
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Fig. 8.11. The ratio of L^, at 850 fim to A/jsm is shown for a sample of low-z spiral 
and starburst galaxies from Dale et al. (2005) and Clements et al. (2010). The 
average and median values for the sample are shown by horizontal lines. 



850 /xm measurements exist with good signal-to-noise and for which global 
fluxes were estimated. 

The 850 fiia fluxes were then converted to specific luminosity -Li^(85o) > using 



= 10-23 47Tdi S^(Jy) = 1.19 X 10^' 5^(Jy) 4(Mpc). 
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^.17) 



Figure 8.11 shows the ratio Li^^^f^/MisM as a function of L,yg^g where 
-^ISM = -^Hi + ^H2- Based on this plot, we then adopt as a working value 



1 X 10^° erg/s/Hz/M0. 



Define this mean value as 



.18) 



asso = 1 X 10^0 erg/s/Hz/M© 



.19) 



Evolution of star formation and gas 



27 



For high-redshift observations, 

observed : {i/obs S^^^J = restframe : ( J^vest^P^ ) , (8.20) 



with 



Therefore, 



-'rest = I^obs X (1 + Z). (8.21) 
Invest Lp . 



(1+^) 4;y^ (8-22) 



-^(i+^Kbs 
with 

= «850 ((1 + z) 33Q°(^'^J ^ MsM. (8.23) 



Then, 

1 +_f (1 + ^)^ob: 

where the 350 GHz is the frequency corresponding to 850 /xm. 
Normalising to Mjsm = 2x 10^° Mq with /3 = 3.8, 



5.,,. = 0850 2 X 10 x^-^(l + .) (^^^) (8.25) 



, (mJy) = 1.67 (1 + z)^-« f ^,7^, (8.26) 

"•^^^ 2x1010^ ' V35OGHZ/ (ii(Gpc) ^ ^ 

at z = 0.3, 1, 2 and 3, ^^(Gpc) = 1.5, 6.6, 15.5 and 25.4 Gpc. 



Figure 8.12| shows the new predicted fluxes as a function of redshift for 



both Band 6 (240 GHz) and Band 7 (347 GHz). 

For reference, looking at the ALMA exposure time calculator, with 
7.5 GHz bandwidth in each polarisation and 10.2 min of integration at Cy- 
cle 1 yields lcj = 0.075mJy at 345 GHz, lo- = 0.042 mJy at 240 GHz and 



0.029 mJy at 100 GHz. The expected fluxes on Fig. [87L2| for z = 2 are ~ 1.7 
and 0.3 mJy respectively for 2 x Mq. Thus Band 7 is optimal since the 
expected flux ratio is ~5:1. Band 3 (100 GHz) is not plotted in the figure 
since its expected flux density is 28 times below that of Band 6. 

To compare with the ability to detect CO, we might use the source 
BX691 from Tacconi et al. (2010) which has a = 7.6 x 10^° and 
=3.5 X 101°Mq at z = 2.19 and CO(3-2) = 0.15 Jy kms'^ For a width 
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Fig. 8.12. The expected ALMA Band 6 (240 GHz) and Band 7 (345 GHz) flux 
densities are shown for Mism — 2 x 10^" Mq. 



of 300 kms~^ , this has an average line flux of 0.5 mJy. If the mass is scaled to 
2 X 10^" Mq, then the average line flux is 0.28 mJy. To get 5a or 0.056 mJy 
sensitivity in a single 300kms~^, requires 5 hours! 

In summary, measurement of the R-J tail of the emission ( and using Equa- 
tion 8.26) thus provides an excellent and fast means of determining dust and 
ISM masses in high-redshift galaxies using ALMA. (The coefficient in Equa- 
tion 8.26| was derived empirically from submm observations and therefore 
may be shghtly different than that obtained from the dust opacity shown in 
Fig. [87}) 
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8.2.7 Effective source size 

For optically thick IR sources, one can estimate the effective size of the 
emitting region from 

/ T \ 1/2 
and scaling to ULIRG luminosities, we obtain 

This effective radius is the overall size of the optically-thick region emitting 
the FIR luminosity. In the event that the emission is optically thin, then the 



estimates from Equation 8.28 are of course lower limits. (In Section 8.2.9 



optically-thick, radiative transfer modelling of the dust emission for a r 



dust density distribution is presented. Figure 8.15 shows the effective radius 



and dust temperature for the emitting region producing the majority of the 



emergent flux - for comparison with Equation 8.28 ) 

For local ULIRGs the typical FIR colour temperatures are ~50 K so the IR 
emission radius is ~500pc for 10^^ Lq. This estimate is similar to the overall 
size of the central concentration in Arp220 (see below), indicating that the 
optically thick assumption is not unreasonable. The most luminous SMGs 
observed at high redshift can have -Lir>1O^^L0; for these sources, the 
emission must come from galactic-scale regions, not just a compact nucleus. 

For the Milky Way the FIR luminosity is ~ 10^° Lq and the mean dust 
temperature ~30-35 K; the effective emitting radius from Equation 8.28 is 



'-^100 pc. However, this emission clearly originates from a large number of 

1/2 

separate clouds, and the mean size of each must be ?^(,loud times smaller. For 
example, if the Galactic emission is assumed to be contributed by ~400 IR- 
luminous GMCs, then the effective size of the IR-dominant region in each 
would be ~5 pc. 



8.2.8 Luminosity and SFR estimates from submm continuum 

Estimating the FIR luminosity (and hence the SFR) from measurements on 
the submm R-J tail is an extremely questionable procedure - this hasn't 
stopped observers from routinely doing it! As noted above the R-J flux 
provides a measure of the dust mass weighed linearly by Tj, but inferring 
a total bolometric luminosity requires knowing where the FIR peaks, and 
the fluxes near the peak. Observations near the SED peak can now be done 
using Herschel PACS (Photodetector Array Camera and Spectrometer) and 
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SPIRE (Spectral and Photometric Imaging REceiver) but many of the SMGs 
are subject to source confusion at the SPIRE resokition. In the absence of 
direct observations at the SED peak, one must assume a dust temperature 
(30-50 K) and optical depth, perhaps based on the submm flux. For many 
of the SMGs, FIR luminosities in the range 10^^~^^ Lq have been estimated, 
implying SFRs of several x 10'^ Mq per yr (typically assuming ~ 30-50 K) 
but such estimates must be extremely uncertain since the derived luminos- 
ity will vary approximately as T^"^. Typical ISM masses of the SMGs de- 
rived from CO measurements or the submm continuum are ~ 1-3 x 10^*^ Mq, 
implying that the ISM will be used up in star formation in an implausibly 
short time of ~ lO'^ yr. 

Blain et al. (2003) attempted to derive empirically a scaling between the 
850 /im flux and Ljr based on the apparent from fitting the SEDs of 
local galaxies. Unfortunately, there is large scatter in the correlation. In 
fact, as we saw earlier, the notion of a single is extremely shaky - both 
because there clearly is a range of temperatures and more importantly, the 
apparent (derived from fitting near the FIR peak) is somewhat degen- 
erate with the dust opacity (which also can cause an exponential fall-off to 
short wavelengths). 

To summarise - it should be clear that one cannot constrain without 
observations near the SED peak; and if one has such observations, it would 
be best to simply use them directly to estimate the luminosity. (Of course, if 
the observed source is at redshift 2; > 5, then the submm fiux measurements 
are in fact probing near the restframe SED peak; they can then provide a 
decent luminosity estimate.) In the next section, we model the optically- 
thick dust sources in order to appreciate some of the systematics and the 
range of uncertainty. 



8.2.9 Modelling optically-thick dust clouds 

For internally-heated FIR sources, it is vital to appreciate that the sources 
have essentially two totally independent parameters: the luminosity L of the 
central heating source and the total mass of dust in the surrounding envelope 
- not much else matters! The character of the spectrum of the source (be it 
young stars or an AGN) makes little difference since the primary photons at 
short wavelength are absorbed in the innermost boundary layer of the dust 
envelope. The radius of this inner boundary is set by the radius at which 
the dust is heated to sublimation (see Section 8.2.2). The overall mass of 
dust of course determines the opacity and therefore the radius at which the 
IR radiation can escape, and thus the 'effective dust temperature' which is 
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Fig. 8.13. The IR SEDs of dust cloud models with 10^^ Lq are shown for increasing 
masses of overlying dust and an r~^ density distribution. The dust heating is 
provided by both the central source and secondary photons from warm dust (see 
text). The peak shifts to longer wavelength as the dust opacity increases and the 
R-J tail rises linearly with dust mass. 

observed. We will see below that the model SEDs of the FIR sources can be 
characterised by a single parameter - the luminosity-to-mass ratio, L/M, 
and thus the problem has, in essence, really just one independent variable 
as long as the source structure is simple (e.g., a single source with radial 
fall-off in density and no clumping). The compactness of the dust cloud is 
parametrised by a radial power law density distribution which can be varied 
but for the discussion below I adopt R~^. 

Using the temperature profiles derived above, I have computed emer- 
gent spectra for a source of central luminosity 10-^^ Lq with overlying dust 
masses ranging from 10^"^ Mq (i.e., total ISM masses ~100 times greater or 
10^"^-^ Mq). These parameters are directed towards ULIRGs and SMGs but 
the results can easily be scaled to lower- or higher-luminosity objects. As 
mentioned above, the critical model characteristic is the luminosity-to-niass 
ratio. For this modelling, the dust distribution is taken to vary as R-^ but 
similar results are found with other reasonable power laws. The inner radius 
is taken at = 1000 K, but this is, of course, not critical since the hot dust 
is covered by the overlying colder dust unless the cloud is optically thin at 
short wavelengths. The outer radius was taken at 2kpc - this also is not 
critical since the dust is cold and optically thin to the secondary radiation 
at the largest radii. 
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Fig. 8.14. The variation in the peak wavelength for and the optical depth 
through the cloud at the peak wavelength are shown for varying L/Mism values 
with a density distribution of dust cx R~^. (Models run for power laws of 0, — 1 
and —2 exhibit similar behaviour as long as the clouds are optically thick.) 



Figure 8.13 shows the emergent specific luminosities {vL^y) for the models 
with different enveloping dust masses but constant overall luminosity. Here 
one clearly sees the effect of varying dust mass and opacity. The clouds 
of lower dust mass have non-exponential short-wavelength SEDs due to the 
lack of high dust extinction on the short-wavelength side of the SED peak, 
and hence the hotter dust in the interior is exposed to our view. This con- 
trasts with the higher-mass clouds which show a peak shifted to relatively 
longer wavelength - due to the fact that dust extinction short of the peak 
precludes photons escaping from the hotter interior radii. Figure 8.13| also 
clearly demonstrates the anticipated correlation (see Section 8.2.6) between 
the flux on the long-wavelength R-J tail and the dust mass. 

To quantify some of the spectral characteristics, Fig. |8.14 shows the shift 
in the wavelength of peak emission ( Apeak i left panel) and the optical depth 
at Apeak to the cloud surface (right panel) for the radial shell contributing 
most to the emergent luminosity. The right panel illustrates what was said 
earlier based on analytics - typically the opacity at the peak will be ~1 for 
a large range of overlying dust masses. Also, since the opacity must be ~1 
at the peak, the wavelength of peak emission will be determined by both 
opacity and dust temperature in realistic models. 

Figure 8.15 shows the radius and dust temperature at which the largest 
fraction of the overall luminosity escapes as a function of the luminosity-to- 
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Fig. 8.15. Left and middle panels show the dust temperature and radius of the shell 
in the model producing the largest fraction of the overall luminosity for varying 
L/MisM values. The right panel shows the optical depth at the peak wavelength 
from this shell to the outer cloud surface - illustrating the fact that the emergent 
luminosity is produced at r < 1. In most cases, the output luminosity is in fact 
spread over a fairly broad range of radii. 



mass ratio. The total luminosity is of course spread over a broad range of 
radii except in the most optically thick models (low L/Mjsm values). 

Figure 8.16 shows the ratio of total IR luminosity (at A = 8-1000 /im) to 
the specific luminosity at A = 850 ^m (i.e., a point on the power-law R-J tail). 
It should be obvious from this figure that there is really no good single value 
for the conversion factor from submm fiux density (Arest ^ 150 fxm) to total 
IR luminosity - these models all had the same total luminosity but varying 
dust masses. Thus it is impossible to reliably estimate the total FIR lumi- 
nosity from R-J flux measurements unless there are additional constraints, 
for example on the luminosity-to-mass ratio or the ISM mass (e.g., from CO, 
or dynamical mass estimates). One simply must have measurements close 
to the IR SED peak. 

There are also some limitations or lessons from this modelling. Most note- 
worthy is the polycyclic aromatic hydrocarbon (PAH) emission which arises 
from transiently heated small dust grains. This small grain component was 
included by calculating the UV-visual radiation energy density at each ra- 
dius and inserting the PAH emission following Draine & Li (2001). These 
short- wavelength features require the emergence of flux from regions where 
the grains are hot. In sources showing such features, the dust obscuration 
must be clumpy, or multiple optically-thin and -thick sources must be con- 
tributing. To illustrate this possibility. Fig. 8.17 shows the effect of reducing 
the dust column in 5% of the surface area to a value 10% of the normal col- 
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Fig. 8.16. The variation in the total IR luminosity to rcstframc 850 fiia flux is shown 
for a range of L/Mism values, illustrating the impossibility of estimating the total 
IR luminosity from a single long- wavelength flux measurement. 

umn density - the NIR/MIR becomes much stronger and the siHcate feature 
appears in absorption. 

The presence or absence of detectable PAH or silicates features depends 
on the covering fraction of the overlying optically-thick dust - their detec- 
tion should not be taken as a reliable indicator of one mode of star formation 
(as done by Elbaz et al. 2011) or alteration of the grain abundances since 
their presence can depend simply on geometry and source non-uniformity. 

8.2.10 Summary 

We have developed a very simple model for the FIR emission which can 
be implemented for fast radiative transfer computations relevant to dust- 
embedded luminosity sources. There are several important conclusions: 

(a) For a central luminosity source, the dust temperature profile with radius can 
be modelled with very simply power laws: Tdust oc r~° "^^ for optically thin 
dust and Tjust oc r~°'^ for optically thick dust. In a realistic source, it will 
be optically thin at the very innermost radius, optically thick at intermediate 
radii and optically thin at the outer radii; the temperature profile is then a 
piecewise fitting together of these two power laws. 
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Fig. 8.17. Models with complete covering (dashed line) and with 5% of the area 
having extinction reduced to 10% and 1% of the complete covering model (solid 
line) are shown to illustrate that a slightly clumpy dust distribution is required 
to model the observed short-wavelength hot dust emission. Transient small grain 
heating is included by calculating the UV-visual radiation energy density at each 
radius and inserting the PAH emission following Draine & Li (2001). 

(b) To model the FIR sources, it is vital to recognise that there are two entirely 
independent parameters: the total luminosity and the mass of dust (and to a 
much lesser extent the radial profile of the dust distribution). Since the tem- 
perature profile scales with L, then most sources may be uniquely characterised 
by the ratio L/M^ust- 

(c) Using these temperature profiles, it is then straightforward and quick to cal- 
culate the emergent SED in the IR as a function of the ratio L/Mjugt. 

(d) At the FIR peak, the emission is often optically thick (based on the steep 
fall-off in the SED to shorter wavelengths) but the optical depth will not be 
very large (i.e., rpeak~2-4) since otherwise the radiant luminosity would not 
be able to escape. 

(e) The effective dust temperature for the emergent luminosity can be determined 
approximately from the wavelength of the peak and then used to estimate an 

effective radius for the source. 

(f) Flux measurements on the R-J tail of the SED are uniquely capable for deter- 
mination of the total dust mass (with a factor of two reliability) since the dust 
is optically thin. And if one can assume a dust-to-gas abundance ratio, such 
measures lead to quick and reliable estimates of a distant galaxy's total ISM 
content. 
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(g) Determination of the total luminosity and hence the SFR requires observations 
near the peak at X^cst ^ 100 /im; it can not be reliably estimated using a single 
long wavelength R-J flux measurement (as is often done for SMGs). 

(h) In optically-thick sources (as judged by the steep fall-ofl^ on the short side of 
the SED peak), the MIR emission at A < 20 /im requires that the dust be 
somewhat clumped with incomplete covering in order to see the hot dust and 
the silicate and PAH features. 



8.3 Star formation in galaxies — two modes 

8.3.1 Quiescent or normal mode of star formation 

As discussed in Section [8. 1[ the star- forming GMCs are self-gravitating with 
internal velocity dispersions implying an effective internal turbulent pressure 
typically 100 times the external diffuse ISM pressures. One's physical intu- 
ition should then suggest that disturbances in the external, diffuse ISM will 
have little influence in general on the rate of formation of stars within GMCs. 
With this in mind, one might expect that normally SFRs will simply depend 
linearly on the overall mass of H2 and/or the internal density of the cloud 
and this logic should hold in other galaxies as long as the CMC properties 
are roughly similar. Galaxy type or the location within a galaxy should have 
little influence on what happens within the self-gravitating clouds. 

Observations of CO emission (tracing the mass of GMCs) and star forma- 
tion tracers (FIR luminosity or Ha) have generally shown a linear propor- 
tionality between SFR and H2 mass distributions. For example, the overall 
Galactic Hil region radial distribution shown in Fig. 8.2 is similar to that 
of H2 as traced by CO. On the scale of individual star- forming clouds, one 
might expect sizable fluctuations in the SFR per unit mass since the overall 
Galactic star formation efficiency (SFE) is low with a characteristic star for- 
mation timescale ~ 10^ yr, i.e., much longer than the CMC internal dynam- 
ical timescale ~ 10^ yr. In addition, the GMCs may have somewhat variable 
mean densities which of course alters their internal dynamical timescale. For 
a sample of individual Galactic GMCs (including ones with and without as- 
sociated Hll regions). Fig. 8.18 (left panel) shows an approximately 'linear' 
correlation between molecular gas mass Myiriai and Ljr or the SFR (Scoville 
& Good 1989) for GMCs with masses 10"^ to over IO^Mq. The scatter off 
this relation is a factor ~4. There is also a clear trend for the most massive 
clouds which all have Hii regions, presumably in the spiral arms, to have 
elevated ratios. 

Given the tendency of the spiral arm GMCs to have somewhat elevated 
luminosity-to-mass ratios, it is probably wisest to use the mean Galactic 
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Fig. 8.18. Left panel: the relation between IR luminosity and H2 mass is shown 
for a large sample of Galactic GMCs with (dots) and without Hii regions (crosses), 
compared with the overall Galactic average obtained from dividing the total Milky 
Way IR luminosity by the total H2 mass (Scovillc & Good 1989). M51 and many 
other nearby spiral galaxies have mean ratios similar to that of the Milky Way. The 
right panel shows a similar plot for a local sample of IR bright galaxies (Sanders 
et al. 1991). The ULIRGs at Lir > lO""^^ Lq show much higher luminosity-to-mass 
ratios, indicating that they are converting gas to stars 10-20 times faster. 

luminosity-to-mass ratio (ALq/Mq) to define a 'normal' mode of star for- 
mation. Such a ratio of Lin/Mf^^c^iiLQ/MQ implies a characteristic star 
formation timescale ~ 10^ yr for conversion of molecular ISM into stars. 
The ULIRGs (at Lir>10^2^q, most of which are gas-rich galaxy merg- 
ers) exhibit a ~ 10-20 times higher Lir/Mhj, hence a 10-20 times shorter 
timescale. 

In summary, there is good basis (empirical from the observations and the- 
oretical intuition from what we know regarding the GMC structures) for an 
approximately linear mode of star formation. From local galaxies such as 
the Milky Way and M51, this is quantified at 2±1 Mq per year per 10^ M© 
of H2, or 

SFR = (MhJMq) • lO-'^ Mq yr-^ (8.29) 

8.3.2 Dynamically-driven starburst mode 

At the same time, there are clearly instances, such as those of the IR-bright 
galaxies, galaxy nuclei and the spirals arms and bars of normal galaxies, 
where the SFE is significantly enhanced to a level such that they can be 
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classified as starbursts. Figure 8.18 (right panel) shows a plot of the IR 
luminosities and molecular gas masses for IR-luminous galaxies (LIRGs, 
Sanders et al. 1991) - all of them have considerably elevated luminosity-to- 
mass ratios compared to the average ratios for the Milky Way and nearby 
spiral galaxies like M51. Both the luminous IR galaxies and the individual 
spiral arm GMCs have exceptionally high rates of massive star formation 
(presumably they are also forming low-mass stars). In such regions, the SFE 
(based on the ratio of Lir/Mhj) may be enhanced by a factor ten compared 
to that in Equation |8.29[ 

In both the spiral arms and the ULIRGs, this elevated SFE appears corre- 
lated with non-circular galactic dynamics and high concentrations of dense 
gas. Galaxy merging leads to dissipative deposition of gas to the centres of 
the merging systems and spiral arm streaming motions cause crowding of the 
GMC galactic orbits. As a result of the non-circular motions and concentra- 
tion of gas clouds, cloud-cloud collisions will be more frequent. If such colli- 
sions occur at relative velocity greater than the GMC internal velocities (few 
km/s), the collisions can compress the internal gas mass and significantly ele- 
vate the SFE (Scoville & Hersh 1979; Tan 2000). This high-SFE mode might 
therefore be referred to as a dynamically-driven starburst where the SFE is 



10-50 times that given in Equation 8.29 - occurring in galaxy mergers or 



localised regions such as spiral arms and nuclear bars of normal galaxy disks. 



8.3.3 Star formation 'laws' 

Over the years a number of prescriptions have been proposed to describe the 
rates of star formation within galaxies. Early on, Scoville & Good (1989) 
and Young & Scoville (1991) advocated that the SFR varied linearly with 
the mass of molecular gas - based on observations of Galactic GMCs and the 
similarity of molecular gas radial distributions in nearby, normal galaxies to 
the radial distributions of SFR tracers (e.g., blue light. Ha and FIR). Ex- 
tensive studies (Kennicutt 1998) which included Hi in the analysis then led 
to the well-known Kennicutt-Schmidt SFR law with SFR oc Shj^jj2 ^'^'^ ^ 
threshold cutoff surface density below which the SFR decreased more steeply. 
In my opinion, the atomic gas (Hi) exhibits almost no correlation with the 
SFR tracers (recently shown clearly by Bigiel et al. 2008 and Schruba et al. 
2011) and is not directly relevant to the star formation process in the inner 
parts of spiral galaxies. (The non-integral, ~1.4, dependence of the SFR in 
the Kennicutt-Schmidt star formation law is likely due to including Hi in the 
analysis and is without a direct physical link to star formation except via the 
formation of molecular clouds and due to the inclusion of starburst galaxy 



Evolution of star formation and gas 



39 




Fig. 8.19. Left panel shows the compilation of observational determinations of the 
SFR and total gas surface densities by Krumholz et al. (2012) together with several 
proposed SFR prescriptions. The right panel shows the ahgnment of these data 
improved when the gas surface densities are normahsed by a free-fall timescale, i.e., 
a scalelength, to provide a fixed 'volumetric' law. 



nuclei. As emphasised above, the GMCs are long-lived and stars form from 
them at low efficiency, so the fact that the Hi is needed to form molecules 
does not justify a physical connection to the star formation process.) 

The most recent compilations have been done by Bigiel et al. (2011) 
and Krumholz et al. (2012). The former advocate a simple linear depen- 
dence on molecular gas surface density; the latter propose a 'volumetric' law 
^* = /h2 ^s'^/ts where the factors entering are the H2 mass fraction, the 
efficiency of star formation per free-fall time, the total gas surface density 
and the free-fall time. Their compilation of observational data is shown in 
Fig. 8.19, together with the various star formation laws they evaluated. Al- 
though their proposed volumetric law provides a good fit to a large range of 
data, the evaluation of rg requires that a scaleheight for the gas be assumed 
and this is not easily derived from the observations so the quality of fit is 
largely dependent on what is assumed. 

At this point, I believe the best approach is simply: 1) assume star forma- 
tion occurs only in the molecular gas (i.e., do not include atomic gas in the 
analysis since it is not physically relevant) and 2) assume two modes of star 
formation: a) a quiescent or normal star formation mode which depends 
linearly on the mass or surface density of H2; and b) a dynamically-driven 
mode (relevant to mergers, spiral arms and bars) which is also linear with 
but has a rate constant up to 20 times that of the first mode. 
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Clearly, for distant galaxies it would be nice to have a more precise star 
formation law with dependence on detailed properties of the clouds (such 
as their internal density and the details of the local galactic velocity field 
streaming and dispersion) but this will not be observationally accessible for 
many many years. Let's keep the star formation prescriptions close to the 
relevant observational capabilities! 



8.3.4 Distinguishing normal star formation and starbursts: 
concentration and timescales 

How does one define a starbursting system (as opposed to a normal star- 
forming galaxy) and what are the observational discriminants between these 
two modes? I think the most physically meaningful definition for a starburst 
is a galaxy or area within a galaxy where the rate of conversion of existing 
ISM into stars is such that the ISM will be consumed in a time significantly 
less than the typical Mism/SFR ~ 1 Gyr timescale for local galaxies like the 
Milky Way. This SFE can be observationally characterised by the simple 
ratio Lir/Mhj. 

Alternatively, one might use the specific SFR per unit stellar mass (SSFR 
= SFR/M*) to identify galaxies where the characteristic timescale (1/SSFR) 
is much shorter than the 'likely' age of the galaxy (which is obviously less 
than the age of the Universe at the galaxy's redshift). At redshift z = 1.5-2.5, 
Rodighiero et al. (2011) argue that there is a tight 'main sequence' of 
star-forming galaxies with fairly constant SSFR at all stellar masses (see 



Fig. 8.29); they then identify starburst objects with SSFR significantly above 
this main sequence. Unfortunately, the SSFRs found for the main sequence 
at z ~ 2 are higher by a factor ~15 than those at low redshift so the discrim- 
inant between the two star formation modes is evolving in time - without 
a clear physical basis, be it a higher gas content, smaller stellar mass at 
early epochs or an increasing efficiency/rate in converting existing gas to 
stars. And, of course, if the definition of the main sequence is evolving in 
time, the fraction of galaxies found to be 'bursting' will depend entirely on 
how far off the main sequence a galaxy must be to be so classified. If most 
galaxies at high redshift were, in fact, bursting, then the 'main-sequence' 
SSFR would be elevated as observed and the dispersion would simply refiect 
the range of burst activity. Clearly, it would be most important now to 
analyse the morphologies of the so-called main-sequence galaxies to classify 
the merging/interacting versus undisturbed percentages. 

Definitions for starburst classification which compare the SSFR to the 
galaxy age or better, using the timescale for exhausting the ISM are more 
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physically based than simply using the dispersion about the main sequence. 
In the near future we can anticipate that ALMA will be measuring gas 
contents for large samples of high-z galaxies. 

Another approach which might be used to identify starburst activity is the 
luminosity density or concentration of the activity. In starbursting systems, 
one expects that a large fraction of the luminosity from the young stars will 
be absorbed by dust and emerge in the FIR and that the luminosity energy 
density within the active region will be high. As discussed in Section [8.2.9 



good IR SED sampling allows one to classify the observed SED according 
to the ratio im/Afdust which is proportional to Lir/Mism- Sources with 
high luminosity density (either AGN or intense starhursts) may also have 
high 'apparent' colour temperatures, providing they are not completely en- 
shrouded by dust which is optically thick out to 100 ^m. One can reasonably 



estimate a characteristic source size using Equation |8.13| and a luminosity 
density from LiR/(47ri?^j^i^i^/3). 



8.3.5 Starbursts in ULIRGs 



In Fig. 8.20 the SEDs, Ljr/Mhj mass ratios and fitted dust temperatures 
are shown for the complete sample of 20 ULIRGs at z < 0.1 from Sanders 
et al. (1988) and Sanders et al. (1991). The SEDs peak at A > 120 /im 
with dust colour temperatures 30-50 K and the FIR-to-H2 mass ratios are 
30-200 Lq/Mq. As noted above, these ratios are a factor 10-20 higher than 
the luminosity-to-mass ratios in local spiral galaxies, indicating a 10-20 
times shorter timescale for converting ISM to young stars. 



8.3.6 Arp220 a prototypical ULIRG 

In the spirit of 'one example can provide more understanding than many 
generalisations', I think it is very worthwhile to look in detail at current 
observations of Arp 220. Of course not all ULIRGs are the same as Arp 220, 
but one can readily see the essence of the starburst mode in this object. 

Arp 220, at 77Mpc, is one of the nearest and the best-known ultra- 
luminous merging system (i^s-iooo/im = 1.5x10^2^0). To power the IR 
output, star formation must be occurring at a rate of ~102 Mq yr~^. Visual 
wavelength images reveal two faint tidal tails, indicating a recent tidal in- 
teraction (Joseph Sz Wright 1985), and high-resolution ground-based radio 
and NIR imaging show a double nucleus (Baan & Haschick 1995; Graham et 
al. 1990). The radio nuclei are separated by 0.98arcsec (Baan &: Haschick 
1995), corresponding to 350 pc. Millimetre- wave imaging provides a unique 
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Fig. 8.20. The SEDs, L^JMh^ mass ratios and fitted dust temperatures are shown 
for the complete sample of 20 ULIRGs at z < 0.1 from Sanders et al. (1988) and 
Sanders et al. (1991). 



capability to probe such starburst nuclei in dusty LIRG/ULIRGs since the 
dust is optically thin at long wavelengths. Virtually all of the ULIRGs ob- 
served so far have massive concentrations of molecular gas in the central few 
kiloparsec (Scoville et al. 1997; Downes & Solomon 1998; Bryant &: Scoville 
1999; Tacconi et al. 1999). 

Arp220 has been imaged at high resolution in the 2.6 mm CO line (Scov- 
ille et al. 1991), 3 mm HCN (Radford et al. 1991), 1.3 mm CO (Scoville et al. 
1997; Downes & Solomon 1998; Tacconi et al. 1999) and most recently in the 
0.9 mm CO and HCO+ lines (Sakamoto et al. 2009). The C0(2-l) line emis- 
sion exhibits two peaks separated by 0.9arcsec, and a larger inclined disk of 
molecular gas (Scoville et al. 1997; Downes & Solomon 1998; Sakamoto et al. 
1999; Tacconi et al. 1999). At 0.5arcsec resolution (Sakamoto et al. 1999), 
the CO and 1.3 mm continuum imaging shown in Fig. |8.2l] reveals counter- 
rotating disks of gas in each of the IR nuclei. The kinematic data require 
very high mass concentrations in both of the nuclei. The nuclear disks are 
counter-rotating, consistent with the notion that the most complete and vio- 
lent merging should be associated with counter-rotating precursor galaxies in 
which there is naturally angular momentum cancellation during the merging. 
The masses of each nucleus are apparently dominated by the molecular gas 
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Fig. 8.21. The merging galactic nuclei of Arp220 are shown in 0.5 arcsec resolution 
imaging of the C0(2-l) and dust continuum emission (Sakamoto et al. 1999). The 
four panels show: a) continuum-subtracted C0(2 1), b) the CO mean velocities, c) 
the 1.3 mm dust continuum, and d) the total CO emission including the extended 
molecular emission. These data resolve the two nuclei (separated by 1 arcsec E-W or 
350 pc) and the crosses indicate the 1.3 mm dust continuum peaks. In the molecular 
gas emission, the two galaxy nuclei have counter-rotating disks as seen in the upper- 
right image. 

- a common finding of tlie ULIRG galaxy studies (Bryant &; Scoville 1999). 

The nucleus of Arp220 has a total molecular gas mass Mh2 ~ 9 X 10^ Mq 
within 300 pc (Scoville et al. 1997; Downcs &: Solomon 1998; Sakamoto et 
al. 1999; Tacconi et al. 1999). This H2 mass within 300 pc is 3-4 times that 
of the entire Milky Way (most of which is in the ring at 4-8 kpc) . 

Although Arp 220 has enormous masses of molecular gas, one should not 
visualise this gas being contained in self-gravitating GMCs like those in 
the Galaxy. The high brightness temperature of the observed CO emission 
(Scoville et al. 1997) requires that the gas distribution be continuous, i.e.. 
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filled nuclear disks, not a cloudy medium. The mass density within the 
disks implies a visual extinction Ay ~ 2000 mag perpendicular to the disks! 
Clearly, studying such regions can only be superficial in the optical/NIR 
and must rely on much longer- wavelength observations. 

The western nucleus in Arp 220 exhibits hard X-ray emission (Iwasawa 
et al. 2005) and very high-velocity CO emission. The inner dust emission 
source seen at A ~ 1 is extremely compact with 10^ Mq within i? < 35 pc 
(Downes &: Eckart 2007). Downes & Eckart argue that the inferred dust 
temperature of ~175 K implies a higher radiation energy density for its heat- 
ing than could be provided by a compact starburst and therefore that this 
nucleus may harbour a luminous and massive black hole. 



8.3.7 An aside: Sgr A* - an extraordinary ISM 

Given the very limited spatial resolution in observations of distant nuclei 
such as Arp 220, it is instructive and cautionary to look at the properties of 
the ISM in our own Galactic nucleus despite the fact that it is not currently 
very active. The massive black hole associated with the radio source Sgr A* 
has a mass of 4 x 10^ Mq - derived from the motions of stars within the 
central 1 pc (Genzel 2006; Chez 2007). A circumnuclear disk (CND) at ra- 
dius out to ~ 3 pc is seen with both ionised gas and dense molecular clumps 
on the exterior. Figure |8.22 shows the ionised gas as probed in the Hi Pa 
line at 1.87 ^m (greyscale; Scoville et al. 2003) in a structure termed the 
mini-spiral and the dense molecular gas as probed in the 3 mm HCN line 
(contours; Christopher et al. 2005). The molecular gas is in a clumpy ring at 
1-3 pc radius. The ionised gas in the mini-spiral may be infalling or outflow- 
ing; the ionised gas at the edge of the disk in the southeastern arm probably 
arises from photoionisation at the inner edge of the molecular cloud ring. 

The total mass of the molecular gas within 3 pc of the black hole is var- 
iously estimated at 1-5 x 10^ Mq and has an orbital time ~ 10^ yr. Given 
that the black hole is only ~ 10 times more massive it is clear that very lit- 
tle of the nearby ISM will actually make it to the black hole, unless we are 
viewing an extremely atypical epoch for Sgr A*. In any case, it is certainly 
the case that the ISM here is not being accreted on an orbital timescale 
since the luminosity of Sgr A* is <C 10^ Lq. 

Just as was the case for Arp 220, our intuition about the properties of the 
ISM in the disk of our galaxy (i.e., self-gravitating GMCs) clearly cannot be 
transferred directly to the ISM in the Milky Way nucleus. The molecular 
gas clumps in the CND have extraordinary properties! Their sizes are <1 pc 
but densities 10^~® cm~^ (Christopher et al. 2005). (By contrast, a typical 
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Fig. 8.22. The circumnuclear disk (CND) in the Galactic nucleus disk in ionised 
gas (greyscale) as imaged in the Hi Pa line at 1.87 /Ltm (Scoville et al. 2003) and 
the dense molecular gas (contours) seen in HON emission (Christopher et al. 2005) . 
The ionised gas is distributed in a 'spiral' pattern and at the inner surfaces of the 
neutral molecular clumps. The densities in the molecular gas are ~ 10''"^ cm~'^(!), 
orbiting at 1-3 pc radius. 



Galactic 10^ Mq GMC has a diameter 40 pc and mean density ~300cm~^, 
see Section 8.1.3). The higher densities are required near SgrA* if the 
clumps are to be stable against tidal disruption, given their proximity to 
the central point mass. If the dust-to-gas ratio in the clumps is similar to 
the standard Galactic value, then the inferred column densities of 10^^ cm~^ 
imply an extinction of Ay ~ 10^ mag through each clump. Since the molec- 
ular ring/torus is clearly not appearing edge-on, its angular momentum is 
not perpendicular to the Galactic plane. The ring therefore probably re- 
sulted from the accretion and disruption of a single cloud at non-zero height 
out of the Galactic plane, rather than smooth accretion over time from the 
larger nuclear disk. If such a cloud arrived at the point where it was tidally 
disrupted with a non-zero height the resulting fragments might end up as 
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Fig. 8.23. A very crude model for the disk in gas-rich merging systems hke Arp220 
is shown. We assume a central point mass (black hole) and a uniform gas surface 
density disk extending to 100 pc radius. The turbulent velocity dispersion within 
the disk is taken to be 50 km/s (based on Arp220; Scoville et al. 1997). The disk 
thickness as a function of radius for hydrostatic equilibrium is shown in the right 
panel and the mean gas density in the left panel (solid red line) . The critical density 
required for stability of a self-gravitating object (e.g., a cloud) in the disk is shown 
by the dashed (blue) curve. The right panel shows the extraordinarily high gas 
densities expected for these very simple assumptions and underscores the fact that 
self-gravitating clouds should not form in the inner disk due to tidal disruption. 



clumps in an orbital plane inclined to the Galactic plane, as is observed for 
the CND. 

At present there is no evidence of star formation within the molecular 
clumps despite their very high densities and extraordinarily short internal 
dynamical timescales ~ 10^ yr. However, the very high extinctions suggested 
above would preclude detection of embedded young stars and the only means 
of detecting them might be via radio detection of ultra-compact Hii regions 
or maser emission. The unusually high densities of the molecular clumps 
must reflect the fact that at lower density they would not be tidally stable 
and therefore could not exist for long (see Section 8.3.8). 



8.3.8 Nuclear starburst disks 

As a first step toward understanding the massive gas disks like that seen in 
Arp220, we consider an extremely simple model with uniform gas surface 
density out to radius 100 pc. At the centre of the disk, we assume a point 





mass of 4 X 10^ Mp) - either a central black hole or a nuclear star cluster 
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Fig. 8.24. The left panel shows the Toomre Q stability parameter (Q < 1 unstable) 
for the disk parameters in the previous figure and the right panel shows the Jeans 
mass for collapse/self-gravitation assuming a sound speed of 50 km/s and neglecting 
tidal shear. 



The disk is assumed to be self-gravitating and in hydrostatic equilibrium 
with a sound speed of 50 km/s, representative of the turbulent velocity dis- 
persion measured in Arp 220. Although this is meant only as an illustrative, 
gedanken, experiment, it clearly shows that the conditions are inescapably 
different from those of star-forming clouds in the Galactic disk. 

In Fig. 8.23, the disk thickness and mean gas densities are shown as a 
function of radius out to 100 pc. The thickness of the disk is typically less 
than 10 pc and the mean densities are everywhere > 10®cm~^, getting up 



to 10^^ cm~'^ near the centre. The right panel of Fig. 8.23 also shows the 
gas density required for tidal stability (dashed line) - clearly demonstrating 
that it will be difficult to form stable clouds within the central 10 pc ra- 
dius. Within this region, we should expect a more continuous gaseous disk 
structure, not a cloudy ISM as found in the centre of our Galaxy. 

The gravitational stability of the disk can also be seen from the Toomre 
Q parameter shown in Fig. 8.24 (left panel). Within the inner 10 pc, Q ^1 
and it is very hard to form self-gravitating structures in the ISM, such as 
clouds. On the other hand, outside 10 pc radius such clouds might form (if 
they are not disrupted by cloud collisions) but their masses are required to 
be very large, exceeding 10^ Mq. 

At the inner boundary of the nuclear disk (i?~ 1-5 pc), one expects that 
the dust will sublimate and any ISM inside this radius will be dust-free and 
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therefore have much lower radiative opacity, permitting radial accretion in- 
wards. On the other hand, at the sublimation boundary, there will be a 
strong radiation pressure gradient in the outward direction, tending to push 
the inner edge of the dust accretion disk to larger radii. The physics of 
this transition zone should indeed be interesting - the material pushed out- 
wards will intercept orbit disk gas with higher specific angular momentum, 
and the dust sublimation interface will also be Rayleigh- Taylor unstable, 
perhaps leading to radial streams of accretion flow. I leave these details to 
the students to work out! 



8.3.9 Maximum-rate starbursts - the dust Eddington limit 

For Arp220, the estimated nuclear SFR of 100-200 Mq per year is at least 
a factor of ten less than the maximal rate, obtained by dividing the ISM by 
the orbital timescale of ~ 2 x 10^ yr. One very effective means of restrict- 
ing the star fromation (or even ejecting the dense ISM) is via the radia- 
tion pressure on dust of the central starburst (or AGN) luminosity. For a 
self-gravitating gas and dust mass, the effective 'Eddington' limit is approx- 
imately 500 Lq/Mq, similar to the overall mass-to- light ratio measured in 
the Arp220 nucleus (IO-'^^Lq/S x 10^ M©). For higher luminosity-to-mass 
ratios, the ISM is blown out of the region. 

The significant role of radiation pressure in high-mass star formation and 
starbursts for feedback through radiation pressure on buildup has only re- 
cently been appreciated. The outward radiation pressure in a starbursting 
cloud core or galactic nucleus will dominate self-gravity at radius R when 

> ^, (8.30) 



A-kR^c R^ 

where (k) is the effective radiative absorption coefficient per unit mass. Al- 
though the original stellar radiation is primarily UV and visible, the dust in 
the cloud core absorbs these photons and re-radiates the luminosity in the 
IR. The 'effective' absorption coefficient takes account of the fact that out- 
side the radius where Ay ~ 1 mag, the luminosity is at longer wavelengths 
where the dust has a reduced absorption efficiency. For the standard ISM 
dust-to-gas ratio (Draine & Li 2007) , Ay = 1 mag corresponds to a column 
density Nh = 2 x 10^^ cm~^ and 

(k)=312-^ cmV\ (8.31) 



where Aeff (i?) is the absorption coefficient- weighted mean wavelength of the 
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radiation field at radius R and for simplicity, we have adopted a A^^ wave- 
length dependence for the absorption efficiency. 



Combining Equations 8.30 and 8.31 , we find that the radiation pressure 



will exceed the gravity when 



For Aeff ~ 3-10 /im, Aog/Ay is approximately 10 and the radiation pressure 
limit is ~5OO-L0/Mq. Thus, for luminosity-to-mass ratios exceeding this 
value radiation pressure will halt further accretion. For a forming OB star 
cluster, this luminosity-to-mass ratio is reached at about the point where 
the upper main sequence is first fully populated, i.e., a cluster with approx- 
imately 2000 Mq distributed between 1 and 120 Mq. A similar 'Eddington' 
limit applies for any dust-embedded starburst region (see Scoville et al. 2001; 
Scoville 2003; Murray et al. 2005; Thompson et al. 2005). It turns out that 
much larger-scale nuclear starbursts such as that in Arp220 have approx- 
imately reached the same empirical limit of 500 Lq j Mq , suggesting that 
nuclear starburst activity may also be regulated by a balance of self-gravity 
and radiation pressure support. 



8.3.10 AGN - starburst: observational connections 

Evidence has also accumulated for an evolutionary link between merg- 
ing ULIRGs and UV/optical quasi-stellar objects (QSOs) as suggested by 
Sanders et al. (1988): similar local space densities for ULIRGs and QSOs; 
continuity of FIR SEDs smoothly transitioning between the two classes; the 
occurrence of AGN-like emission lines (Veilleux et al. 1999) and significant 
point-like nuclei (less than 0.2 arcsec - Scoville et al. 2000) in 30-40% of the 
ULIRGs; and the association of both ULIRGs and some QSOs (MacKenty &: 
Stockton 1984; Bahcall et al. 1997) with galactic interactions (see the review 
article by Sanders & Mirabel 1996). Whether the entire QSO population 
had precursor ULIRGs (implying that galactic merging is the predominant 
formation mechanism for AGN) is certainly not yet settled and at lower 
AGN luminosities, much of the activity is probably not directly associated 
with galaxy merging or interactions. 

Two possible scenarios linking the ULIRG and luminous AGN phenomena 
are: 1) that the abundant ISM, deposited in galactic nuclei by merging, fuels 
both the nuclear starburst and feeds the central black hole accretion disk; or 
alternatively, 2) the post-starburst stellar population evolves rapidly with 
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a high rate of mass return to the ISM in the galactic nucleus - leading to 
sustained fuelling of the black hole (e.g., Norman Sz Scoville 1988). 

8.3.11 AGN - starburst: theoretical connections 

The largest potential sources of fuel for AGN are the nuclear ISM (until it 

has been cleared out) and the mass loss associated with normal stellar evo- 
lution of stars in the galactic nucleus. Most recent discussion has assumed 
the former link; so I would like to briefly discuss the latter possibility here 
- partially as a stimulus for further investigation by students. For a nuclear 
starburst population, approximately 20% of the initial stellar mass will be 
lost via red giant mass- loss winds within the first 2 x 10*^ years (Norman 
&: Scoville 1988). (The mass return associated with supernovae is probably 
significantly less.) Scoville Sz Norman (1988) and Scoville & Norman (1995) 
examined the fate of this mass-loss material - specifically to account for both 
the broad emission lines (BELs) and the broad absorption lines (BALs) seen 
in AGN. Since the stars will not be disrupted by AGN feedback once the 
nucleus becomes active, the stars can provide a long-term supply of fuel to 
power the AGN. 

The physics of the dust shed by the stars in their stellar winds, par- 
ticularly its evaporation, is critical to determining whether the mass-loss 
material accretes inwards to an accretion disk or is blown outwards by ra- 
diation pressure. This consideration leads naturally to a division of the 
central cluster environment into: 1) an inner zone (r < Ipc) where the dust 
is heated to above the sublimation temperature and the lower opacity of the 
dust-free gas allows it to fall inwards, and 2) an outer zone (r > 1 pc) where 
the dust (and gas) survives and is driven outwards at velocities up to 0.1c 
by radiation pressure. 

It is interesting to note that the density required for tidal stability of a 
BEL cloud at 1 pc radius from the black hole requires that the clouds must 
have extremely high internal densities. At 1 pc from a 10^ Mq black hole, 
the Roche limit density is 3 x 10^^ cm~^. This provides a strong argument 
for the stellar (rather than interstellar) scenario to account for the BELs, 
using the stellar masses to bind the emission line gas. 

8.4 Evolution of galaxies at high redshift 

Over the last 15 years extensive surveys of high- redshift galaxy evolution 
(e.g., GOODS, UDF, COSMOS, AEGIS and CANDELS) with space- and 
ground-based telescopes have dramatically increased our understanding of 




Fig. 8.25. A compilation oi z — 4 and 6 UV LFs for galaxies (Capak et al. 2011). 
The curves are Schechter function fits by Bouwens et al. (2007) for z = 4 (blue) and 
6 (black) galaxies and the points are derived for Lyman alpha emitters (LAEs) and 
Lyman break galaxies (LBGs). 

galaxy evolution in the early universe. All of these go very deep (typi- 
cally > 27-30 mag AB) but cover very different size areas - most also have 
excellent multi-wavelength, ancillary data. The largest survey (COSMOS) 
has over a million galaxies with photometry and photometric redshifts at 
z = 0.2-6 and the deepest (UDF) now has detections at z = 6-8, probing the 
first 1 Gyr of cosmic time. Major evolution is seen in the galaxy properties: 
stellar mass, luminosity, size and SFR. 



8.4-1 Luminosity and Mass Functions 

Figure |8.25| shows a compilation of recent determinations of the z = 4-6 
luminosity functions (LFs; Capak et al. 2011). Evolution of the UV LF 
(probing the distribution of star-forming galaxies) is clearly seen with the 
number densities increasing at all luminosities as one goes to lower redshift 
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Fig. 8.26. The stellar mass function for galaxies at z = 0.1-4 from Marcliesini et al. 
(2009). 



and there is apparent steepening of the low-L power law going to higher 
redshift, i.e., more low- luminosity galaxies contributing. If this persists to 
higher redshift then it can be argued that the reionisation of the universe 
at z ~ 10 must have been produced by relatively low-luminosity galaxies 
(Robertson et al. 2010). At z>6, the samples are very small, ~ 20-50 ob- 
jects with virtually no spectroscopic confirmation. Even at z = 6, one can 
see large scatter in the different determinations of the LF, shown as points in 



Fig. 8.25 It is noteworthy that the high-L portion of the LF is very poorly 



constrained (see Fig. 8.25) and could even be a power law rather than the 
Schechter function (exponential fall-off) used in fitting. 

At lower redshifts, comparatively good samples exist for the stellar mass 
functions of galaxies (see Fig. 8.26). Strong growth in the number density of 
galaxies at all masses is seen but the largest increase at late epochs (low z) 
occurs in the lower- mass galaxies (see, e.g., Ilbert et al. 2010). This 'down- 
sizing' of evolution at later epochs is also seen in the galaxies with star 
formation and AGN activity. 
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Fig. 8.27. Large-scale structures mapped by the projected density of galaxies in 
redshift slices (Az = 0.02-0.2) in the COSMOS survey field (Scoville et al. 2007a; 
Scoville et al. 2012). Over 200 significantly overdense regions are seen over the 
redshift range z = 0.1-2.55. 

In principle, the buildup of stellar mass in galaxies as a function of time 
should be consistent with the evolution of star formation activity within 
galaxies in each mass bin. Discrepancies between these two might indicate 
the rate of galactic merging (assuming the masses and SFRs are reliable). 
This has been explored by Drory & Alvarez (2008) who find reasonable 
consistency between the stellar mass buildup and the SFR evolution - pro- 
viding the SFRs drop first in the highest-mass galaxies and galaxies typically 
undergo ~ 1 major merger after z = 1.5 (see also Reddy 2011). 

8.4-2 Environmental correlations 

In ACDM simulations for the early Universe, it is the most massive, highly 
biased, structures which form earliest. Such structures will become the lo- 
cations of massive galaxy clusters with the most massive galaxies built there 



54 



Nick Z. Scoville 



first. The COSMOS survey (Scoville et al. 2007b) was specifically designed 
to probe a sufficiently large area on the sky (2 square degrees corresponding 
to > 50-100 comoving Mpc at z> 0.5) that the full range of environmental 
densities could be sampled at all redshifts. This enables both the map- 
ping of the large-scale structure and the investigation of the correlation of 
galaxy evolution with environment. High- accuracy photometric redshifts 
and a large sample of spectroscopic redshifts enable the separation of galaxy 
structures along the line of sight. Figure 8.27 shows the 200 overdense re- 
gions of galaxies seen in COSMOS at z = 0.1-2.55 (Scoville et al. 2012). 

Using the environmental densities shown in Fig. |8.27[ I show the percent- 
age of early-type galaxies (with SED corresponding to E-Sa galaxies) as a 
function of density and redshift in Fig. 8.28 (left panel). As is well known 
from many studies, the early-type galaxy fraction increases systematically 
to a fraction ~ 50% at z = 0. However, Fig. 8.28 (left panel) also clearly 
isolates the differential evolution associated with large-scale structure den- 
sity - the early-type galaxies form first in the highest-density regions of the 
large-scale structure. 

The evolution of the SFR density (SFRD = SFR per comoving volume) 
rises a factor 20 from z = to a peak at z = 2-3. An important question 
is: in which environments is the star formation occurring at each epoch? 
Figure 8.28 (right panel) shows that the SFRD systematically shifts from 
dense to less dense regions of the large-scale structure as time progresses. 
Peng et al. (2010), using both COSMOS and SDSS data, show that the 
quenching of star formation activity in galaxies has two separable terms: 
one dependent on galaxy mass and the other on environmental density. 

Recently, Rodighiero et al. (2011) have attempted to quantify the frac- 
tion of star formation occurring in starbursts versus normal-rate star for- 
mation at z = 1.5-2.5. They make use of a BZK colour-selected sample of 
galaxies with low SFR and samples from GOODS and COSMOS, probing 
with Herschel and Spitzer the moderate- and high-FIR luminosity galaxies 



see Fig. 8.29). The BZK galaxies are used to define a 'main sequence' of 



star-forming galaxies at this epoch. The outliers above the main sequence 
(mostly IR-detected) with more than four times the SSFR are then taken to 
be galaxies undergoing burst-mode star formation. Rodighiero et al. (2011) 
suggest that only 2% of the galaxies by number are undergoing starbursts 
and this population contributes ~ 10% of the total SFRD at this epoch. 
This is an interesting approach to the issue, but it should be clear that the 
adopted definition of the main sequence and its spread will entirely deter- 
mine the starburst percentages. In addition, some of the galaxies within 
the main sequence could well be merger-starbursts of originally lower-SSFR 
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Fig. 8.28. Left panel: the early-type (E-Sa) galaxy fraction is cle arly correlated 
with both redshift and large-scale structure density from Fig. 8.27 (Scoville et al. 
2007a; Scoville et al. 2012). The early-type galaxies form first in the highest- 
density environments. Right panel: the percentage of the overall SFRD is shown as 
a function of redshift and large-scale structure density. The dominant environments 
for star formation shift to lower density at later epochs (Scoville et al. 2012). 



galaxies. As noted earlier, it would be advantageous to have a physically- 
based definition of bursting rather than simply picking outliers as bursts. 
To take an extreme example, the latter definition would clearly be wrong if 
essentially a large fraction of the population were, in fact, bursting. 



8.5 Modelling star formation at high redshift: 
same modes but different frequency 

It is now well established that the activity associated with both star forma- 
tion and AGN increases dramatically (20 x) out to a peak at z~2-3 using 
SFRs determined by the UV continuum, FIR and radio continuum and that 
beyond z = 3 there may be a gradual decline in the SFR density, although 
the measurements become increasingly difficult at these redshifts. 

At high redshift, one might expect that the same two modes pertain, yet 
their relative importance could be quite altered. On the one hand, high- 
redshift galaxies should have higher gas mass fractions than the typical 
5-10% of low-redshift spiral galaxies like the Milky Way - leading to higher 
rates of star formation associated with the quiescent linear mode of star for- 
mation. But at the same time, one expects a greatly elevated rate of galactic 
merging/interaction, increasing the frequency of the starburst mode. Which 
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Fig. 8.29. The SFRs of three samples of galaxies at z — 1.5-2.5 and stellar masses 
are shown (Rodighiero et al. 2011). The samples include BZi^T-selected galaxies 
(black dots), GOODS (cyan) data and COSMOS (red) galaxies with Hershel PACS 
FIR detection. They use the BZK galaxies to define a 'main sequence' of star- 
forming galaxies at this epoch and then classify the outlier (mostly IR-detected) 
galaxies at > 4 times the SSFR on the main sequence as starburst-mode galaxies. 

of these dominates is not at all clear without a numerical simulation to track 
their relative importance, keeping track of the rare and brief mergers. 

8.5.1 Cosmic evolution: 
M^, and Mjsm (^nd star formation luminosities 

In order to test the framework developed above against the observed cos- 
mic evolution of galaxies, I developed a simple Monte Carlo simulation - 
including: star formation in the quiescent and merger-driven burst modes. 
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merging of dark matter haloes and their contained galaxies, and accretion 
of fresh gas from the external large-scale structure environment. The sim- 
ulation starts at z = 6 with a population of 10 million dark matter haloes 
having a mass function approximating that seen in simulations for z = 6 
(e.g., Heitmann et al. 2010). Each halo is started at 2; = 6 with a galaxy 
of baryonic mass equal to the dark matter mass times the universal baryon 
fraction, with 90% of the baryons being gaseous ISM and 10% stellar mass. 
I then let the galaxy and dark matter halo population evolve with 50Myr 
timesteps down to z = 1. At each timestep, gas is converted to stellar mass 



at the quiescent SFR given by Equation 8.29 At each timestep, haloes 



are also randomly sampled for having undergone a merger with a proba- 
bility weighted as (1 + ^)^-^M^^{25 (pakhouri & Ma 2010). This merger 
rate is normalised such that 2% of the haloes merge per 50Myr at z = 6. 
The secondary galaxy for each merger was also selected randomly from 
the galaxy population weighted by (mass ratio)^^'^ (Fakhouri &: Ma 2010). 
For those galaxies selected to merge, the SFE was increased by a factor of 
10-50 (but only for one timestep). Gas accretion to the galaxy halo was 
taken as M = 6.6(Mhaio/10^2)^-^^(l + zf-'^^ x 0.165 (Dekel et al. 2009) for 
Mhaio<lO^^M0. For lar ger-mass haloes, we assumed simply that the ac- 
cretion was cut off - either by the accretion shocks or by AGN feedback - 
the former is a departure from Dekel et al. but some reduction of accretion 
is required in order to have the massive objects become gas-poor ellipticals 
at modest redshifts as shown observationally. 

Minor non-critical details which were included were that: 1) the effective 
accumulation of stellar mass was taken to be 70% of the integrated star for- 
mation (i.e., assume 30% of the stellar mass is recycled eventually in mass 
loss), 2) the accretion of external gas to star- forming galaxy was delayed by 
1 Gyr after it accreted to the halo boundary (to account for the infall time) 
and 3) during starbursts, ISM mass was shed from the galaxy at a rate equal 
to the SFR (only for the burst mode). The star formation luminosity (moti- 
vated to model the IR luminosity function) was taken very crudely as the to- 
tal luminosity from stars formed in the last 50 Myr plus that from young stars 
in earlier timesteps reduced by a factor of two in each timestep. Specifically, 
the luminosity associated with recent star formation was taken to be 10^° Lq 
per Mq per year of star formation, based on observations of local galaxies. 



8.5.2 Need ISM replenishment by accretion 

This simplistic model was remarkably useful to explore critical aspects of the 
evolutionary scheme involving quiescent and burst-mode star formation with 
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Fig. 8.30. The mass and luminosity functions for galaxies at z = 2.5 (left) and 
z = 1 (right) derived from a Monte Carlo simulation including galaxy merging 
and associated starburst activity but no gas replenishment via accretion from 
the external environment. The dashed line is the original (z = 6) galaxy mass 
distribution (gas+stars) . The original gas supply is exhausted far too quickly, even 
by z ~ 2. 



a reasonable gas accretion hypothesis. Figure 8.30 shows the evolved mass 
functions of stars and ISM gas and the star formation luminosity at z = 2.5 
and 1. For this figure, the simulation included merging galaxies and their as- 
sociated starbursts but was without accretion of gas from the environment. 
In this case, the ISM runs down at a rate given by local Universe star forma- 
tion laws and the original gas content of the haloes is exhausted by z = 2.5 to 
a level much less than that seen either at z = 2 or in present-epoch galaxies. 



8.5.3 ULIRG starbursts account for high-L tail 



The simulation shown in Fig. 8.31 includes gas accretion as formulated above 
but without the starburst activity associated with galaxy merging. Now, 
the ISM and stellar mass functions exhibit characteristics similar to those 
observed at z = 2.b and 1, i.e., < 50% gas mass and 10%, respectively. How- 
ever, even at z = 2.5, the star formation luminosity function is lacking the 
power-law tail at the high-luminosity end. To reproduce the power-law tail, 
the starburst activity associated with galaxy merging is needed, as shown 
in Fig. 8.32 (The low-mass and low-luminosity ends of the distribution 



functions rise more steeply than is observed since no effort was made to 
model the behaviour there - this shallow slope is often attributed to star- 
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Fig. 8.31. Similar to Fig. |8.30| except that gas accretion at halo masses below 
10^^ M0 is included, following the prescription of (Dekel et al. 2009). In this 
simulation, the starburst activity associated with merging has been removed and 
the high end of the luminosity function naturally must follow the exponential form 
of the mass distribution. 
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Fig. 8.32. Here both accretion and starburst activity during galaxy merging are 
included, giving a reasonable qualitative match to the observed high-end mass and 
luminosity functions. 



formation winds driving galactic mass-loss at velocities above the escape 
velocity of lower-mass galaxies.) 
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8.6 Conclusions 

In these lectures, I have attempted to summarise a physically intuitive ap- 
proach toward understanding the characteristics of star formation in high- 
redshift galaxies making use of our intuition gained from the low-2; Universe. 
I have devoted considerable effort towards developing an understanding of 
the molecular clouds (GMCs) in which most low-redshift star formation oc- 
curs in normal galaxies. 

I have also provided a summary of the physics of IR emission from opti- 
cally thick dust clouds. Since there is very little systematic treatment from a 
theoretical point of view of the observed IR SEDs, I spent considerable effort 
to develop a systematic approach. Most of these IR sources are centrally 
heated by a starburst or AGN, so the theoretical framework for the dust 
heating and radiative transfer can be greatly simplified and distilled. Inter- 
preting the observed SEDs is particularly prone to mistakenly evaluating a 
dust temperature from the wavelength of the SED peak when in reality a sig- 
nificant range of dust temperatures is contributing to the emergent flux and 
the peak is equally determined by the overall dust opacity (since obviously 
the radiation at a given wavelength cannot escape if the opacity is too high 
at that wavelength). I underscore also the point that the long- wavelength 
R-J tail fluxes should be used to derive dust (and ISM) masses - not total 
IR luminosities (as is commonly done for submm galaxies). 

Given the known internal structure of the star-forming GMCs, it is rea- 
sonable to expect that star formation activity will have two modes: 1) a 
quiescent mode in which stars form in GMCs at a rate roughly proportional 
to the mass of H2 gas (accounting for the linear correlations observed be- 
tween the CO emission and star formation tracers such as FIR over the disks 
of local galaxies including our own) and 2) a dynamically-triggered starburst 
mode in which the rate per unit gas mass is elevated by factors of 10-50 ac- 
counting for the preferential formation of Hii regions in spiral arms, despite 
the presence of GMCs throughout the disk, and the high luminosity-to-gas 
ratios seen in local ULIRGs. For the starburst mode, collisions of GMCs 
may signiflcantly increase the internal structure/density of the H2 gas and 
stimulated star formation due to supernovae and expanding Hll regions. 

At high redshift, the relative importance of these modes will be shifted, al- 
beit in different directions, by two changes: 1) higher gas-mass fractions (in- 
creasing the quiescent mode) and 2) higher rates of galactic merging (increas- 
ing the starburst mode). Some observational determinations of the galaxy 
merger rates at high redshifts have shown very discrepant results for the evo- 
lution. However, my recent analysis of evolution in the frequency of galaxy 
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close pairs from COSMOS data suggests merger evolution as (1 + z)^'^, in 
good agreement with the dark matter halo merger rates seen in simulations. 

At high redshift the conditions in the gas may be quite different from 
those in the star- forming clouds at low redshift. It has been speculated 
that the stellar IMF might be top-heavy in starbursts or low-metallicity en- 
vironments. This is not an insurmountable problem - observations of the 
4000 A break may be used to constrain the mass of low-mass stars. Until 
such observations arc done, my conservative opinion would be to avoid the 
'last refuge of scoundrels' since very little hard evidence is found in the local 
Universe to support such variations. 

To test whether the ingredients discussed above provide a reasonable ba- 
sis for understanding high-redshift galaxy evolution, I show the results of a 
Monte Carlo simulation starting at z = 6 with 10 million haloes distributed 
in mass and with merging rates as found in ACDM simulations. Starting 
with galaxies for which the gas content is 90% of the baryonic mass, the 
galaxies are evolved including both star formation modes and gas accretion. 
The simulations clearly show the need for gas replenishment through accre- 
tion from the external environment; otherwise, the z = 2 and present-epoch 
gas contents are far too low. 

This need for accretion is simply a reflection at high-z of the well-known re- 
quirement of gas replenishment for the Milky Way. Specifically, for the Milky 
Way the SFR ~ 3 Mq per year and the present gas content is ~ 3 x 10^ Mq, 
implying an exhaustion timescale of only 1 Gyr. In the case of local galaxies, 
it does not appear that the accreting gas is Hi since the high velocity Hi 
clouds do not constitute a sufficient infiux. More likely the inflow is in the 
form of diffuse ionised hydrogen (Hii) for which imaging with high sensitivity 
to low surface brightness emission is required. Clearly, this is an important 
direction for future observations. 

Another conclusion from the simulations is that the starburst mode trig- 
gered by galactic merging is necessary to account for the high-luminosity 
power-law tail of the IR luminosity functions at high redshift. The quies- 
cent mode cannot do this since the luminosity generated would simply reflect 
a scaled version of the galaxy mass function which is exponentially falling 
at the high mass end. 
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